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ABSTRACT 

Wolf-Rayet stars (WRs) are evolved massive stars, and the relative number 
of WC-type and WN-type WRs should vary with metallicity, providing a sen- 
sitive test of stellar evolutionary theory. The observed WC/WN ratio is much 
higher than that predicted by theory in some galaxies but this could be due to 
observational incompleteness for WN-types, which have weaker lines. Previous 
studies of M33's WR content show a galactocentric gradient in the relative num- 
bers of WCs and WNs, but only small regions have been surveyed with sufficient 
sensitivity to detect all of the WNs. Here we present a sensitive survey for WRs 
covering all of M33, finding 55 new WRs, mostly of WN type. Our spectroscopy 
also improves the spectral types of many previously known WRs, establishing 
in one case that the star is actually a background quasar. The total number of 
spectroscopically confirmed WRs in M33 is 206, a number we argue is complete 
to ~5%, with most WRs residing in OB associations, although ~2% are truly 
isolated. The WC/WN ratio in the central regions (<2 kpc) of M33 is much 
higher than that predicted by the current Geneva evolutionary models, while the 
WC/WN ratios in the outer regions are in good accord, as are the values in the 
SMC and LMC. The WC/WN ratio and the WC subtype distribution both argue 
that the oxygen abundance gradient in M33 is significantly larger than found by 
some recent studies, but are consistent with the two-component model proposed 
by Magrini et al. 
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Group — stars: evolution — supergiants — stars: Wolf-Rayet 

1. Introduction 

Wolf-Rayet (WR) stars are the final evolutionary stage of high mass stars before they 
erupt into spectacular Type Ibc supernovae (Meynet et al. 2011). While in this stage, stars 
can be thought of as "bare stellar cores" characterized by highly evolved surface elements 
(Meynet & Maeder 2005) and high mass-loss rates of ~1O~^M0 year~^ (Abbot & Conti 1987, 
Nugis et al. 1998, Crowther 2007). These strong stellar winds create the broad emission lines 
that characterize a WR star's spectrum (Cassinelli & Hartmann 1975, Conti & de Loore 
1979). Depending on the emission lines present in the spectrum, a WR star can be one of 
two types: WN (nitrogen rich) or WC (carbon rich). Both types contain large amounts of 
helium with little to no hydrogen while WN stars contain nitrogen and WC stars contain 
both carbon and oxygen (Massey 2003 and references therein). 

According to modern stellar-evolutionary theory, the most massive stars (and thus, the 
stars with the highest mass-loss rates) will first lose enough of the outer layers to reveal the 
H-burning products, helium and nitrogen, (WN stage) before continuing to remove enough 
material to reveal the He-burning products, carbon and oxygen (WC stage). However, less 
massive stars will only reach the WN stage. This is known as the Conti Scenario (Conti 
1976, Maeder & Conti 1994). It follows that WC stars will be more common relative to 
WN stars in galaxies with high metallicitie£], as the higher metallicities will result in higher 
mass-loss rates. The WC phase will also last for a larger fraction of the He-burning phase 
(see, e.g., Meynet & Maeder 2005). 

Current evolutionary models do a good job at predicting the relative number of WCs 
to WNs at low metallicities (SMC, LMC), but seem to fail at high metallicities. This 
is well illustrated in Figure 11 of Meynet & Maeder (2005) for the Geneva evolutionary 
models, where the data for the WC/WN ratios comes from Massey (2003) and references 
therein. Various efforts have been made to solve this problem, i.e., Eldridge & Vink (2006) 
and Meynet et al. (2011). However, the fault may not be that of the models but of the 
observations themselves. Conti & Massey (1989) showed that the strongest optical line in 



^Throughout this paper we use "metalhcity" to reahy mean the oxygen abundance. This shorthand is 
justified in studying the evolution of massive stars, as oxygen is one of the primary drivers of the stellar 
wind. See discussion in Massey et al. (2003). Fortunately, too, the oxygen abundance is (in principle) 
straightforward to measure in extragalactic HII regions. 
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WCs (C III A4650) has a median line strength a factor of 4 higher than for the strongest 
optical line in WNs (He II A4686), and thus WNs are harder to detect. So, while a survey 
may be complete for WCs, it may not delve deep enough to detect all of the WNs. This 
would skew the WC to WN ratio and thus make the evolutionary models appear incorrect. 
At the same time, surveys of a few small regions within a galaxy may skew the results to 
either higher or lower WC/WN values as the population may not be well mixed in age. With 
this in mind, a galaxy-wide complete survey must be undertaken to measure the WC to WN 
ratio before the models can criticized and improved if needed. 

In this paper, we present the first survey of the entire disk of M33 with sufficient 
sensitivity to detect even the weakest-lined WNs. M33 is particularly interesting to examine 
because of its uncertain metallicity. Studies of WR stars in other galaxies have shown that 
the WC to WN ratio is quite sensitive to metallicity. For example, compare the SMC 
(log(0/H) + 12 = 8.13) and the LMC (log(0/H) + 12 = 8.37) (Russell & Dopita 1990). 
Massey et al. (2003) found the SMC ratio to be 0.09 while Breysacher et al. (1999) found 
the LMC ratio to be 0.22. Previous WR studies of M33, such as Massey & Conti (1983) and 
Massey & Johnson (1998) show a decrease in the ratio of WCs to WNs as the galactocentric 
distance increases (center ratio of 0.6, outer ratio of 0.2). Additionally, Massey & Johnson 
(1998) found a change in the WC spectral properties; WCs in the outer regions showed broad 
lines while WCs in the inner regions showed narrow lines. These findings are consistent 
with a M33 having a significant metallicity gradient, such as that found by Garnett et al. 
(1997). However, recent studies have found that there is little or no metallicity gradient 
based upon HII regions with measured [OIII] A4363 electron temperatures (Crockett et al. 
2006, Rosolowsky & Simon 2008, Magrini et al. 2007, 2010, Bresohn 2011). Nevertheless, 
Magrini et al. (2007, 2010) argue that the metallicity in the central region of M33 is much 
higher based on the abundances in young massive stars. Thus, by presenting a complete WR 
survey of M33, we hope to further investigate the metallicity of M33 as well as compare our 
results to evolutionary theory. Subsequent studies will revisit the issue of the WC/WN ratio 
of the low-met allicity starburst galaxy ICIO, and present the first galaxy-wide deep survey 
for WRs in the relatively metal-rich disk of M31. 

In Section [2] we describe how we identified candidate WR stars through the use of 
photometry and image subtraction techniques. In Section [3] we discuss the spectroscopic 
confirmation process. In Section H] we present an updated list of all of the currently known 
WR stars in M33. In Section |5] we evaluate our completeness and discuss our results and 
finally, in Section E] we summarize our findings and list future goals. 
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2. Identification of Candidates 

We have begun a new survey for Wolf-Rayet stars among the star-forming galaxies of 
the Local Group (beyond the Magellanic Clouds), using the basic filter system introduced by 
Armandroff & Massey (1985). The goal of the present survey is to remove the bias towards 
the detection of WC stars discussed extensively by Massey & Holmes (2002). 

Owing to their strong emission lines, Wolf-Rayet stars are relatively easy to detect. Al- 
though slitless spectroscopy (objective prism) surveys were effective at finding Wolf-Rayet 
stars in the Magellanic Clouds (Azzopardi & Breysacher 1979a, 1979b), crowding and faint- 
ness (compared to the night sky) compromises their use for distant galaxies (e.g., Bohannan 
et al. 1985). However, the use of narrow-band interference filter imaging has proven highly 
effective. This technique was first used by Wray & Corso (1972) who surveyed two 8' di- 
ameter fields near the center of M33 using photography through an interference filter that 
included both the C III A4650 line and the He II A4686 emission lines. The same fields were 
photographed through a continuum filter, and 25 WR candidates were identified by visually 
blinking the on-band and off-band exposures, most of which were subsequently confirmed 
as WRs, predominately of WC type (see Table 9 of Massey & Johnson 1998 and references 
there in). Armandroff & Massey (1985) described a three-filter set that was designed to be 
more sensitive to finding WRs, and to distinguish between WNs and WCs: one filter ("WC") 
is centered on the C HI A4650 line, another ( "WN" ) is centered on He II A4686 (strong in 
both WNs and WCs), and a continuum filter ("CT") is centered at A4750. Each filter has a 
band- width of approximately 50 A. Armandroff & Massey (1985) used CCD imaging through 
this set to identify WR candidates in NGC 6822 and IC 1613, as well as two small fields in 
M33 that overlapped with those of Wray & Corso (1972). Spectroscopy confirmed that the 
Wray & Corso (1972) imaging had missed many WNs, not unexpected given the difference 
in fine strengths between CHI A4650 in WCs and He II A4686 in WNs. Massey & Johnson 
(1998) used a version of this set to image eight 5^2 x 5^2 fields in M33, and Massey & Duffy 
(2001) imagined the SMC through such a filter set and identified some additional WNs that 
had been missed by objective prism surveys. 

For the work here on M33 (and our survey of other Local Group galaxies that will be 
described in subsequent papers) we used the Kitt Peak 4-meter Mayall telescope and the 
Mosaic CCD camera to image through a similar set of filters. The data on M33 were obtained 
on 2000 Sept 21, 23, and 24 (and allowed to properly age). Three overlapping 36' x 36' fields 
were used to cover a 0.32 deg^ area, extending ~8.5 kpc in radius from the center and 
covering all of the optical disk. Each field was imagined using three 300s exposures in each 
filter, with the telescope dithered slightly between the exposures to fill in the gaps between 
CCDs in the Mosaic camera. The seeing varied from I'.'l on the center and southern fields 
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to I'.'S on the northern field. 

Flat-fielding of the images was achieved by exposures of an illuminated uniform screen 
attached to the dome, and bias frames were obtained every night. Reductions followed the 
same procedure as for the Local Group Galaxies Survey (LGGS) data described by Massey 
et al. (2006), and included corrections for "cross-talk" between pairs of chips that shared the 
same electronics. The only additional step was that since these filters were non-standard, 
higher-order spatial distortion terms had to be derived before transforming to a linear plate 
scale. In analyzing the images, we treated each of the three dithers as separate exposures; 
i.e., we did not attempt to "stack" the images as this would compromise their photometric 
integrity as emphasized by Massey et al. (2006). 

As we now describe, we used two techniques with these images in order to identify 
stars that were "significantly" brighter on the on-band rather than the continuum. The first 
involved photometry of individual stars while the second was based upon image subtraction 
techniques. 

2.1. Photometry 

Crowded field photometry was used to identify candidate WR stars by comparing the 
magnitude differences of WC - CT and WN - CT to the magnitude uncertainties based on 
poisson statistics and CCD read noise. Instrumental magnitudes of every detected star were 
obtained using IRAF's point-spread function (PSF) fitting routine "DAOPHOT" (Stetson et 
al. 1990) and the magnitude differences were then computed. Next, the number of measured 
stars and their positions in clusters were matched between the WN, WC and CT frames. 
Stars were added or subtracted by eye to obtain matches between the frames, and the PSF- 
fitting routines were re-run. This allowed for the detection of stars in even the most crowded 
regions of M33. Stars with magnitude differences larger than -0.10 mag and a significance 
level greater than 3a were considered candidates. 

For the photometry candidates, there were 1693 WR candidates with counterparts in 
the LGGS found to a low significance level (3(t). If one were to use the absolute magnitude 
of M33 to scale from the LMC {My = —18.9 and My = —18.5, respectively; van den Bergh 
2000), one might expect to find ~190 WR stars in M33, or about 30% more than the 141 
WRs Massey & Johnson (1998) list. Thus, while we do expect to find a few more WR stars, 
the number of discoveries should be on the order of tens, not hundreds. Therefore, this 
large number of candidates exceeds the expected number of WRs by about a factor of 100. 
This gave us some indication that we needed to be a little more selective with our list of 
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candidates and that the instrumental errors were probably underestimated. Since we knew 
that the Massey & Johnson (1998) survey was deep enough to detect even Of-type stars, 
we instead considered the known M33 WRs. On average, they were detected at a level of 
25(7, and only six (out of 137 unambiguously identified) were found at a significance level 
less than 10a. The lowest significance one of these was found at 5.2a, while the next least 
significance known WR was found at the 7.1a level. Thus we decided to keep only those 
photometric candidates with significance levels > 5a; there were 382 such objects including 
many of the known ones. 

2.2. Image Subtraction 

We complemented the photometry with the use of image subtraction, a technique that 
has achieve much success in the detection of supernovae (see, for example. Yuan et al. 2010). 
Since WR stars will be brighter in either the WN or WC filters than in the CT images, 
simply subtracting the CT image from the WN image should reveal an image displaying 
only the WN stars. Similarly, subtracting the CT image from the WC image should reveal 
the WC stars. However, seeing variability often turns this simple idea into a much more 
complex problem. To solve this, we used the Astronomical Image Subtraction by Cross- 
Convolution program which employs point-spread fitting techniques and cross- convolution 
to help overcome these barriers (Yuan & Akerlof 2008). 

After running the image subtraction program on each of the three dithered images in 
each field, we were able to proceed with identifying potential WR candidates. We first took 
the median of the images. This reduced the anomalies created by bad pixels as well as 
erroneous results due to poor image subtraction when dealing with bright foreground stars. 
By visually blinking through the images and comparing the stars on the WN or WC image 
and the result image to the original CT image, we were able to remove many false positives. 
As Figure [1] shows, the WR candidates we identified were of one of two types: either it 
was easy to visually identify the WR star since the magnitude differences were large, or 
it was virtually impossible to distinguish any magnitude difference between the two filters. 
Still, even stars with seemingly negligible magnitude differences showed up prominently as 
potential WR stars in the resulting image. 

A study of all three M33 fields using image subtraction yielded 330 WR candidates, 225 
(68%) of which matched with candidates identified by photometry. Additionally, 129 of the 
152 (85%) known WR were identified using this method. Why the remaining 23 confirmed 
WR stars were not identified is discussed in Section 15.11 
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3. Spectroscopy 

In order to determine which of our WR candidates were real, we used the 6.5-m MMT 
with the 300 fiber- fed spectrometer Hectospec (Fabricant et al. 2005). The instrument 
covers a 1 deg diameter field of view, well matched to our survey area, and we prepared 
two configuration files with different field centers to cover the maximum number of new 
candidates. The spectra were taken on UT 2010 Oct 9 and 11 (program NOAO 2010B- 
0260) with the 270 line mm~^ grating blazed at 5000 A, and covered the wavelength range 
3700-9000 A. Beyond ~7500 A the spectra may be contaminated by second-order blue light, 
but this has no impact on our project, which relied upon data in the 4000-6000 A range for 
identifying and classifying WR stars (i.e., from N IV A4058 through C IV A5812). The fibers 
have a diameter of 250 jj^m (1'.'5 on the sky), and resulted in a spectral resolution of 6 A (5 
pixels). Wavelength calibration was by means of a He-Ne-Ar comparison arc taken in the 
afternoon, with a wavelength shift (in pixel space) apphed from the prominent A8399 OH 
night-sky line. Flat-fielding was accomplished by an exposure of an illumination calibration 
screen, with a through-put correction for each fiber determined by its location within the 
field. Complete details of the reduction procedure can be found in Drout et al. (2009). Fibers 
in each configuration were assigned to carefully selected blank sky positions. 

We made use of the UCAC2 (Zacharias et al. 2004) to select suitable guide stars known 
to have negligible proper motions. This catalog is on the global International Celestial 

Reference System (ICRS), while the LGGS coordinates are tied to the USNO-Bl.O catalog 
of Monet et al. (2003), a predecessor to the ICRS. Using the UCAC2 catalog we applied 
the following transformations to the LGGS coordinates to make them match the ICRS: 

aucac2 = Q!LGGS " 0!030 and 5ucac2 = ^LGGS - 0'.'14. 

In assigning priorities for the two field configurations, we concentrated first on the 
objects found both by the image-subtraction method and photometrically but were not 
previously known. The second priority targets were those found with a high significance 
level (> 6a) from the photometry comparison. Third priority were those stars found by 
image subtraction but not found by the photometry. Fourth priority were lower significance 
photometry candidates not found via the image subtraction. And the last priority were 
previously know WRs, as we did not want to waste fibers. (Not all of the known M33 WRs 
were included.) As we shall see, obtaining high quality spectra for some of the "known" 
WRs proved illuminating. 
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3.1. "Losers" and Non-WR Star "Winners" 

In our first glance at ttie spectra, we were gratified to find a large number of newly 
discovered WR stars. But even a cursory examination revealed that many of our WR 
candidates were in fact late-type stars. The reason was obvious from the spectra: there 
is a TiO molecular feature at A4760 within the CT bandpass that caused these stars to 
appear brighter in the emission line filters than in the continuum. This feature only becomes 
pronounced at spectral types of M3 and later (see Fig. 14.1 of Jaschek & Jaschek 1990), and 
the effect was somehow overlooked when the second author originally designed the filter set, 
probably as there were only a few stars this late in the Jacoby et al. (1984) stellar library 
(see Armandroff & Massey 1985). The problem equally affected the image subtraction and 
photometric candidates. Adding to our chagrin, of course, was the fact that such stars could 
have been readily removed from the sample simply by consulting the colors listed in the 
LGGS. Indeed, all of the stars for which we see TiO have B — V > 1.1, with the exception 
of one star with unrealistic colors (i.e., B — V ^ —1.3) due to crowding issues. 

We list in Table [1] the total number of stars in each category and number observed, 
(more to the point) the number of non-red stars. We see that the vast majority of the stars 
detected with both methods (image subtraction and photometry) turn out to have emission, 
which we consider to then be a "winner" , even if the star is not actually a WR. The other 
9 stars (plus the 44 found only by image subtraction) invariably are bright stars, usually 
bright supergiants. These stars, although interesting in their own right, will be discussed 
elsewhere as part of our survey of the less evolved massive star population. Of the non- 
emission photometric candidates, the losers covered the entire range of possibilities, from 
0-type dwarfs (with He II A4686 in absorption) to mid F-type foreground dwarfs. The most 
likely explanation are mismatches in the number of stars identified on the continuum and 
emission-line images, despite our best efforts. If a tight grouping was photometered as 5 
stars on the CT image but only as 3 or 4 stars on the WN or WC images, then there will be 
spurious high significance candidates. 

Not all of our "winners" were WRs or Of- type stars (Table [2]). In fact our spec- 
troscopy discovered two previously unrecognized quasars (QSOs), one of which had been 
previously classified as a Wolf-Rayet star! The spectra of the QSOs are shown in Figure |2j 
J013445.02-f304928.0 has a redshift of z = 1.99, while J013322.09+301651.4 had a redshift 
of z = 2.84. The latter was classified as a WN star by Massey & Conti (1983) based upon a 
much poorer signal-to-noise spectra, which showed only a single broad feature near A4686; 
see their Figure 2. The line is actually Lyman eg Note the rich Lyman „ absorption fo. 



'To a carpenter, everything looks like a nail." 
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est short-wards of this hne evident in our spectrum. Both objects are x-ray sources in the 
XMM- Newton Serendipitous Source Catalog and also appear in the Atlas of Radio/X-ray 
Associations by Flesch (2010). Perhaps some day these background objects can be used 
to probe the ISM of M33, although their faintness {V = 19.4 and 20.1) would make such 
high dispersion studies difficult at present^. We also rediscovered the "hot LBV" candidate 
JO 13459.47+303701. 9 found by Massey et al. (2007b). Several of its relatively weak (-1 to 
-2 A equivalent width) Fell and [Fell] lines were in the WC and WN filters, but the CT 
bandpass happened to be free. One of our objects (J013331. 23+303334. 2) is an HII region 
with strong He II A4686 emission, apparently not previously known. The object was not 
considered an He II A4686 emission nebula by Kehrig et al. (2011) as it was too stellar (M. 
S. Oey, private communication 2010). 

Perhaps the most interesting two non-WR "winners" are shown in Figure El The stars 
J013423.00+304536.5 and J013351.46+304057.0 both show strong lines of H and He I, along 
with weak [Nil] AA6548, 6584 lines but no sign of the [OIII] emission that would characterize 
a normal HII region. JO 13423. 00+304536. 5 shows weak lines of He II. The lines are all 
blue-shifted since the object has a velocity of -230 km s~^. Given the star's location within 
the plane of M33, this velocity is about as expected (Zaritsky et al. 1989; Drout et al., in 
prep). The emission- line spectrum resembles that of some symbiotic stars (Munari & Zwitter 
2002) although the late-type giant absorption features one expects in such objects in the far 
red are not present (Kenyon 1986). (The impression of absorption in the figure are due to 
slight over-subtraction of the night-sky spectrum.) The star is fairly bright {V = 20.47, 
B — V = 0.27, U — B = —1.06, all according to the LGGS) and if we adopt an average 
reddening of E{B — V) = 0.12, then the absolute visual magnitude My is —4.5, which is also 
consistent with a symbiotic. J013351. 46+304057.0, on the other hand, shows very broad, 
weak emission that extends from 4600 A to 4720 A, and appears to include some He II 
A4686 as well as He I A4713. This broad feature is why the star made it into our list, but the 
spectrum is not that of a WR star. Its spectrum also includes the nebular [S II] AA6716, 6731 
(located in the figure just to the right of He I A6678). The fiux ratio of [SII] A6716 to A6731 
is about 1.4, requiring that these lines originate in a low density region (Osterbrock 1974). 
It is quite a bit brighter than that of J013423.00+304536.5, with V = 17.73, B -V = 0.08, 
and U — V = —0.93, according to the LGGS), and hence My = —7.2, which is too bright 
to be a symbiotic star. We leave the nature of both objects unresolved, a subject for future 
researchers. 



■^http://wnimssc-www.star.le.ac.uk/Catalogue 

''We are indebted to David James for this suggestion. 
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3.2. Newly Found and Previously Known Of-type and WR Stars 

In Table [3] we list the Of-type and WRs that we observed spectroscopically, including 
previously known WRs. We measured the equivalent widths (EWs) and full widths at half 
maxima (FWHMs) of the "discovery" lines, He II A4686 for the WNs, and the C III A4650 
and He II A4686 blend in WCs. We also measure the EWs and FWHMs of C HI A4650, and 
C IV AA5801,5812 feature (hereafter referred to as C IV A5806) if present. For WCs, the 
very strong C HI A4650 feature is blended with (and dominates over) He II A4686. 

Classification of these stars was, for the most part, straightforward. The primary classi- 
fication scheme is derived from Smith (1968a, 1968b) for WN3-WN8 and WC5-WC9 stars, 
with the extension to WN2 and WN9 and WC4 by van der Hucht et al. (1981). See Conti 
(1988, p. 13) for a complete summary. In earlier investigations of the WRs stars located 
beyond the Magellanic Clouds, the signal to noise of the spectra was often too poor to make 
more than a crude division into "early" or "late" types (Massey et al. 1987). The signal- 
to-noise we achieve here is modest by some standards (typically 40 per spectral resolution 
element), but the result is that our spectra often resemble textbook examples of their class. 
We illustrate the spectra of a few representative examples in Figures |4] and |5] covering the 
range of subtypes we observed. 

The higher signal-to-noise afforded by the 6.5-m aperture of the MMT and good through- 
put of Hectospec allowed us to detect the weak absorption lines in stars that might otherwise 
have been classified as a weak-hned late-type WN star. For instance, J013333. 58+304219. 2 
had been called a WN-type Wolf-Rayet star by Massey & Johnson (1998), where it was listed 
as C22, but our superior spectrum here shows that it is actually an 04 If star, with both He I 
A4471 and He II A4200,4542 absorption. There is nebulosity present (as evident from [OIII] 
AA 4959, 5007 and the [Nil] lines flanking Ha), and the Balmer lines show a narrow emission 
component, further complicating the interpretation. We show the spectrum in Figure |6l as 
we believe this is the earliest supergiant yet found in M33 (i.e., based upon an inspection of 
Table 17 of Massey et al. 2007b and Table 9 of Massey et al. 2006). Although an equivalent 
width of He II A4686 of -10 A is often taken as the dividing line between Of-type stars and 
WNs, there are stars whose EWs fall in the grey area, and J013333. 58-1-304219. 2 is one such 
star. Abbott et al. (2004) has similarly reclassified J013344.68+304436.7 (OB66-25) as an 
Of-type star rather than the WN8 classification offered by Massey et al. (1996), and we have 
reclassified J013432.68+304655.4 (N604-WR10 from Drissen et al. 2008) as an Of-type star 
rather than WN6 since we can see absorption. We list the stars that are no longer considered 
WRs in Table H 

The Ofpe/WN9 class (sometimes known as "slash stars") was introduced by Walborn 
(1977) to describe four bright emission-line stars in the LMC with spectral properties in- 
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termediate between Of-type supergiants and WNs. Spectral and photometric variability 
has linked these to the Luminous Blue Variables (LBVs); see Bohannan & Walborn (1989). 
The spectroscopic survey of UV-bright luminous stars in M33 by Massey et al. (1996) dis- 
covered a number of such objects. Using UV-spectroscopy, optical data and model fitting, 
Bianchi et al. (2004) showed that these M33 objects tend to have chemical abundances 
intermediate between O supergiants and WNs. Our survey here detected several of these 
objects, some previously known. These include the peculiar object found by Romano (1978) 
[J013509. 73+304157.3] and classified by Massey et al. (2007b) as Ofpe/WN9, as well as 
J013416.07+303642.1 and J013418.37+303837.0, newly found here. We illustrate the spec- 
tra of these three stars in Figure [3 Note the strong P Cygni features at He I A4471 and the 
broad sea of N II emission blended with He II A4686. 

We also re-observed the Ofpe/WN9 star J013245.41+303858.3 (UIT 008), and were 
astonished to find that its characteristic P Cygni He I A4471 line has changed to a pure 
absorption profile. In Figure [8] we compare our 2010 spectrum with that obtained on 19 
Nov 1995 and published in Massey et al. (1996). Both spectra were obtained with multi- 
fiber instruments, and hence the sky subtraction is not local; thus, nebulosity is not easily 
removed. The spectra are clearly contaminated by nebulosity, as evidenced by the strong 
[OIII] AA 4959, 5007 lines in our 2010 spectrum. (The 1995 spectrum did not extend beyond 
4900 A.) The fiber diameters used for our new observations have half the diameter (1'.'5 vs. 
3'.'0) than those used in 1995, and indeed the Balmer lines are much stronger in the earlier 
spectrum. While it is undecided whether these are primarily nebular in origin, there is no 
question that the He I A4471 profile in 1995 was P Cygni, not absorption contaminated 
by emission. Thus we might have naively classified this star as an 08If and attributed 
the remaining Balmer emission to nebulosity. We retain the Ofpe/WN9 designation here, 
although without He I A4471 in P Cygni this classification is somewhat a matter of curtes}|f|. 

Another remarkable spectrum is shown in Figured The star J013510. 27+304522 would 
normally be classified as a WN2, as its spectrum shows only lines of He II, but no N HI, 
N IV, or N V. What makes this spectrum so peculiar is the presence of extremely strong O VI 
A3811. Similar WC stars are known, and such objects have given rise to the (seldom-used) 
'WO" classification (Barlow & Hummer 1982). Nugis & Niedzielski (1990) argue that this 
line is weakly present in some WNs. But, it is clear from Figure M that the O VI line is 
exceptionally strong in this star. Oxygen is enriched as a result of He burning in massive 
stars, and so we would naively expect to see strong carbon lines, while only very weak C IV 
A5806 is present. Such a star would offer a strong challenge to single-star evolutionary models 



^Note that Crowther et al. (1995) prefer to call such stars WNIO or WNll rather than maintaining the 
original "slash" designation. 
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if it signals high oxygen abundance. 

4. The Current Census of Wolf-Rayet Stars in M33 

We hst in Table S] all of the 206 currently known Wolf-Rayet stars including the 55 new 
WR stars discovered in this survey. We have removed the 3 stars from the 141 "confirmed" 
WRs of Massey & Johnson (1998) that we now consider to be non-WRs (Table |2]), and added 
in the 4 newly found WRs by Massey et al. (2007b), the 7 by Drissen et al. (2008) and the 2 
by Kehrig et al. (2011). Although Massey & Johnson (1998) attempted to introduce a WR 
numbering scheme for M33, we adopted a coordinate based system here, as we were pleased 
to find that all of the known WRs have a counterpart in the LGGS survey of M33, even in 
the crowded HII regions studied by Drissen et al. (2008) with HST. We acknowledge that 
the designations are cumbersome, but at least it should remove any ambiguity about which 
star is which. We include cross references to the Massey & Johnson (1998) numbers, as well 
as some previous IDs, in Table |H Additionally, we give the designations from Drissen et al. 
(2008) even though in many cases their new designations referred to previously known WRs. 

Included in the table are what we consider to be the best spectral types; for instance, 
Abbott et al. (2004) obtained better spectra for a number of the Massey & Conti (1983) 
stars and thus we've adopted Abbott et al.'s spectral types, unless our current spectroscopy 
improves upon it. 

There are admittedly some arguable inclusions in the list. For instance, J013332. 64+304127. 2 
currently shows a spectrum resembling that of a hot LBV star, but we have included it as it 
was spectroscopically WN for many years, as discussed in Massey et al. (2006). Yet we do 
not include J013416.07+303642.1, which Crowther et al. (1997) consider to be a 'WNllh", 
as discussed below. 

We include in the table the value p (p of 1 = 7.53 kpc after adopting a distance of 
840 kpc to M33, following Freedman et al. 1991), the de-projected distance from the center 
of M33 within its plane. These values differ substantially from the older ones we used, 
primarily as the result of using the de Vaucouleurs et al. (1976) Second Reference Catalogue 
/i25 isophotoal radius of 30 f8 rather than the often-quoted Holmberg radius of 25' as the 
normalization factor. We also adopted more modern values for the inclination (56°) and 
position angle (23°), following Zaritsky et al. (1989), although these are quite similar to 
what we had previously used, following Kwitter & AUer (1981). The center of M33 is taken 
to be a2ooo=l''33'"50!89 and (52ooo=30°39'36'.'8. 
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4.1. Cross-identifications 



In creating a single list containing all of the known Wolf-Rayet stars in M33 we were 
occasionally stymied by some of the old identifications. For instance, Massey & Conti (1983) 
identify two WN stars (MC8 and MC9) near the HII region MAI but provided only "ap- 
proximate" coordinates taken from instrument header. Massey & Johnson (1998) retained 
these two stars as "MC8" and "MC9" and repeat the approximate coordinates. In Table 9 of 
Massey et al. (2006), both (!) J013303.19+301124.2 anti J013302.55+301135.8 are identified 
as MC 9 while J013302.94+301122.8 is identified as MC 8. Compounding this, Kehrig et 
al. (2011) list both MC 8 and MC 9 as J013303. 19+301124.2, although in point of fact the 
correct two stars are identified in their Figure 3 (as WR14 and WR15, respectively )|§. To 
resolve this, we successfully resurrected the original spectra of MC 8 and MC 9 taken by 
Massey & Conti (1983). The authors had taken spectra of the two very bright knots in the 
MA 1 HII region, and MC 8 corresponds to J013302.94+301122.8 and MC 9 corresponds 
to J013303. 19+301124.2. (The second star identified as MC9 in Table 9 of Massey et al. 
2006, J013302.55+301135.8, is not a WR star.) Neither MC8 nor MC9 were identified in our 
current imaging study as WRs as the objects were too bright {V = 16.1 and V = 17.2) for 
the image subtraction to work reliably, and the region too crowded for the photometry to be 
believable. Given that the absolute magnitudes of early WNs tends to be —3 or —4, while 
these objects have My of —9 and —8, we can safely conclude that these WRs are actually 
marginally unresolved clusters of early-type stars, perhaps scaled down versions of R136. 
See further discussion in Kehrig et al. (2011). 

Another difficult case was that of J013416. 07+303642.1. This star was "discovered" by 
Massey et al. (2007b) to be an Ha emission line star, and spectrally resembles that of P 
Cygni. It was included in Massey et al. (2007b) as an LBV-candidate, cross-referenced to 
H 108 of Corral (1996). This star turns out to be none other than Humphreys & Sandage 
(1980) B517, the spectrum of which was classified by Crowther et al. (1997) as a Wolf-Rayet 
star of type "WNllh". We would not ourselves consider P Cygni to be a WR star, nor do 
we consider this object to be one, and therefore we do not include it in our table her^. 

Drissen et al. (2008) used a combination of interference filter imaging with HST and 
ground-based spectroscopy to identify WR stars (and candidates) in the giant HII regions 
of M33, i.e., NGC 592, NGC 595, and NGC 604. The superior resolution of HST helps in 



^We are indebted to Carol Kehrig for useful correspondence on the matter. 

^The star had also been observed earlier as Esteban (1994), who had misidentified it as coincident with 
the WC-type WR MC70. Crowther et al. (1997) conclude that the two sources are not the same, and here 
we confirm that MC 70 is J013416. 28+303646. 4, located slightly a few arcseconds east and 4" north of B517. 
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such crowded fields (particularly NGC 604), but by carefully comparing their finding charts 
to the LGGS images, we found that all but two of their confirmed and candidate WRs were 
well resolved on our ground-based LGGS images and appear as unique objects in the LGGS 
catalog. The exceptions are J013333. 75+304134.0 (MC 31) which their HST imaging re- 
solves into N595-WR2A and -WR2B, both of which may be WRs, and J013433.62+304704.6 
(MC 76) which they call N604-WR6, which has an unresolved companion in our LGGS pho- 
tometry. They list cross identifications of two of their stars each for MC 74 (CM 11) and 
MC 75 (CM 12), but this must reflect their uncertainty in the cross-identifications rather 
than evidence of blending given the reasonable spatial separation between objects. We do 
not include their N604-WR2 or N604-WR3 stars in Table [5] as these have not yet been con- 
firmed spectroscopically, and one of their WR candidates did prove to be an 06.5 laf star 
(N604-WR13). In addition, our spectrum of J013432.68+304655.4 (their NGC604-WR10) 
clearly shows He II A4200 and A4542 absorption lines, leading us to call this star an Olf 
rather than a WN. (The difference is presumably due to our higher signal to noise.) Note 
that while their ifS'T-derived coordinates were very precise they are not very accurate, as 
inconsistencies in the HST guide star reference system from one region to another results in 
errors as great as 1'.'5 compared to the ICRS now used in astrometric catalog^. 

Other cross-identification issues were extensively discussed in the appendix of Massey 
& Johnson (1998). 



4.2. OB Association Membership 

Massey & Johnson (1998) list the membership in the OB associations identified by 
Humphreys & Sandage (1980). Here we have extended this to the newly found WRs, using 
their charts to assign the membership indicated in Table |5l Although most of the WRs 
in the tables are in or near OB associations, 42 of the 206 (20%) are not0. Most of these 
stars are nevertheless located along spiral arms, primarily near the center of M33 where the 
photographic work of Humphreys & Sandage (1980) was compromised by saturation. But 
in a few cases (we estimate 2%), these stars are truly isolated. Examples of truly isolated 
WRs are listed in Table [61 At the opposite extreme, we find small OB associations that are 
chock-full of WRs. A good example of this is the OB association 102, which has a diameter 



^Recall that we needed to transform our own LGGS coordinates slightly to place them on the ICRS; i.e., 
Section [3l 

^We include as members the stars in the MAl (Mayall & Aller 1942), even though this is not cataloged 
in Humphreys & Sandage (1980) as it is outside their survey area. 
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about 0'.5 (about 120 pc) but contains 3 WRs. 

4.3. Absolute Visual Magnitudes 

It may be of interest to consider the absolute visual magnitudes (My) of M33 WRs. The 
values are poorly known for Milky Way WRs (due to uncertain distances and reddenings), 
and for the Magellanic Clouds, due to the smaller sample size and limited range of spectral 
subtypes. In Table [7] we list the average My for each spectral subtype, where we have 
removed stars of uncertain type or whose spectrum shows an absorption-line signature. We 
compare these values to those for the Magellanic Clouds. We have treated all of the data 
similarly, except that for the SMC nearly all of the stars have an absorption line signature 
and we have left these in. 

There are three interesting things that we note from these data. First, the standard 
deviations (a) are extremely large for all subtypes that have a reasonable number of stars 
[N) in the sample. The data here are for broad-band V, and hence can include emission 
lines, but similarly large scatter is seen in the narrow-band absolute magnitudes given in 
Table 3-2 of Conti (1988, p. 122). We believe that this reflects the fact that stars of a variety 
of masses and luminosities pass through each subtype; i.e., that WN3 stars (for instance) 
do not come from one mass progenitor, while WN5 stars come from some other. Second, 
it is interesting to note that the My's of the M33 WRs are no higher or lower, in general, 
than that of the same subtypes in the LMC. For galaxies beyond the Milky Way, we expect 
that line-of-sight companions will increasingly contaminate the sample due to crowding, but 
we expect the presence of such companions bright enough to affect the My significantly will 
manifest itself spectroscopically through the presence of absorption lines. The fact that the 
My's are similar for LMC and M33 WRs without obvious absorption signatures confirms 
this. Third, there are 3 WN2s in our "clean" sample, and these have a considerably brighter 
My than the (one) LMC WN2 star. Based upon very small number statistics one might 
otherwise infer that as a subclass these stars are particularly faint visually (Conti 1988, p. 
122), but this may not be true in general. 

4.4. Additional Candidates 

There were only 7 "rank 1" WR candidates that we didn't have a chance to observe 
spectroscopically; we expect 5 of these to actually turn out to be WR stars rather than 
Of-type stars, or (say) quasars. We list these remaining stars in Table [Hi 
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5. Results 

5.1. Completeness 

To understand the completeness of our survey, we must examine any biases created 
during the candidate selection or spectroscopic confirmation processes. Luckily, we are able 
to ignore any bias caused by differential reddening in M33 because the range in E(B-V) is 
modest (0.09 - 0.33, Massey et al. 1995). However, there are a few other aspects we must 
examine. 

While the photometry method of candidate selection was able to identify all of the 
known WR stars, the image subtraction method only identified 129 (or 85%) of the 153 
previously known stars. However, this 15% difference was not due to any shortcomings in 
image subtraction per se - instead, it was due to the Mosaic CCD images themselves. When 
identifying candidate WR stars using the image subtraction method, we required that the 
same star be seen on each of the three dithered images. Thus, only the area of M33 that was 
covered by three images was surveyed. This area works out to be ~85% of the area covered by 
photometry. All 23 known WR stars (except one which we just plain missed) that were not 
identified by the image subtraction technique were located in the "gaps" caused by dithering 
the Mosaic CCD images. Therefore, we can safely assume that while image subtraction did a 
wonderful job of finding new WR stars, we missed finding ~15% of the yet-to-be-found WR 
stars. However, considering we found 55 new WR stars, and this method was complemented 
by the photometry, at the high-end this works out to only be 8 stars or ~5% of the WR 
stars in M33. 

Crowding is often an issue when observing young massive stars because they haven't 
had time to drift away from their birthplaces in OB associations and tightly packed clusters 
(Massey & Conti 1983). By using crowded field photometry, as described in Section [2?T| we 
were able to partially circumvent this issue. However, even though a great deal of effort was 
spent on this, in some cases the process wasn't as successful as in others. Image subtraction 
appears to have done a better job identifying WR stars in clusters since the missed stars 
mostly fell in the Mosaic CCD gaps and not in dense clusters. Thus, the intersection of the 
two resulting sets forms our best list of candidates. 

The most important question we must ask in terms of our completeness is: is our survey 
sensitive to the weakest-lined WNs? As can be inferred from Table |3l the identification 
lines in the weakest-lined WNs are more than 5 x weaker than the lines in the weakest-lined 
WCs. Thus, we will concern ourselves with our WN-completeness and ignore the WCs. It 
has long been known that the line strength of He 11 A4686 of the most extreme Of-type 
stars is comparable to the that found in weakest WN stars (Conti & Frost 1977). As our 
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spectroscopy confirmed, we detected 14 Of-type stars due to their He II A4686 emission, and 
so we might conclude from this alone that our survey was indeed sensitive enough to detect 
even the weakest lined of WR stars. However, Massey & Johnson (1998) argued that in a 
survey such as ours, what is primarily important are the emission-line fluxes. It is clear that 
for two stars of equal brightness, but with different He II A4686 EWs, that it will be harder 
to detect the star with the weaker line. However, it is also true that it will be easier to 
detect the brighter of two stars with the same EWs of He II A 4686, all other things being 
equal. Massey & Johnson (1998) used the V-band flux of a star (based upon its magnitude) 
times the EW to approximate the line flux, as absolute spectrophotometry was not available 
for their data, just as it is not for ours. (Fiber data is notoriously hard to fiux-calibrate 
well.) In Figure [TO] we show the data from the current study. Among the M33 WNs (open 
circles) we see the same trend found by Massey & Johnson (1998): the fainter stars also 
have the weaker line fluxes (i.e., that the EWs are very similar). The solid points denote 
the Of-type stars, and they indeed have a weaker line-fiux for He II A4686 than WRs of the 
same brightness, i.e., they have weaker equivalent widths. We include in this diagram line 
fluxes of the WNs in the LMC (blue x's, from Conti & Massey 1989) as well as for the SMC 
(red asterisks, from Massey et al. 2003). We find that our new survey has detected WNs 
with lines as weak or weaker than those of WNs in the Magellanic Clouds, galaxies whose 
WR content is essentially completely known (see discussion in Massey 2003). 

It is also important to point out the quality of our spectroscopy, which is better even 
than that of Massey & Johnson (1998). With our spectra, we can easily distinguish between 
Of-type stars and WN+absorption (WN stars whose emission is diluted by the continuum of a 
companion). Massey & Johnson (1998) were unable to do this reliably because their signal to 
noise was not alwasy able to detect the absorption. As an example, see J013419. 58+303801. 5, 
whose equivalent width is less than 10 A. Generally, we use -10 A as the division between 
an Of-type and a WN and so at first pass we may have not classified this star as a WN. 
However, the signal to noise in our new spectrum is 75 per 6 A spectral resolution element, 
and with that we can detect quite a few other He II lines in emission, thus making this star 
a WN as such emission is not found in Of-type stars. With such high signal to noise, we 
are confident that if we identified a WR candidate and took spectra of it, we were able to 
correctly classify the star. This was not always possible in previous studies. 



5.2. The Metallicity of M33 



As discussed in Section [T|, the relative number of WCs and WNs (WC/WN) is quite 
sensitive to metallicity. Massey & Conti (1983) and Massey & Johnson (1998) found that 
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in M33 this ratio decreased with galactocentric distance consistent with what was beheved 
to be a strong metalhcity gradient at the time (Kwitter & AUer 1981, Zaritsky et ah 1989, 
Garnett et al. 1997). Recent studies have found a much shallower slope for the oxygen 
abundance (Crockett et al. 2006, Magrini et al. 2007, Rosolowsky & Simon 2008, Bresolin 
2011) based upon HII regions with [OIII] A4363 measurements, although Magrini et al. 
(2010) have emphasized that this selection may significantly bias the sample towards low 
metallicity, particularly in the inner regions. 

Figure [TT] shows the distribution of the 206 known WC (red +'s) and WN (blue x's) stars 
throughout M33. It is clear that, even with our more complete data, there is a significant 
difference in the relative numbers of these stars, with the WCs concentrated towards the 
center. The green ovals represent distances of p = 0.25 (1.9 kpc) and p = 0.50 (3.8 kpc) 
within the plane of M33. In Table M we list the number of WCs and WNs within the 
corresponding regions. The stated errors on the ratios are based upon the assumption of 5% 
incompleteness. We plot these WC/WN values against p in Figure [T2l 

How well do these values agree with the measured metallicities for M33? In Fig- 
ure [T^left) we show the medley of recent determinations of the oxygen abundance plot- 
ted against galactocentric distance. (For easy comparison we have normalized the distance 
within the plane of M33 in the same way as for the WR data, i.e., p = 1 corresponds to 
7.53 kpc) First, we see that there is good agreement for p > 0.4 (3 kpc). Inwards of this 
region, studies have been usually been compromised by the scant number of HII regions 
that have been analyzed (Crockett et al. 2006) or large scatter (Rosolowsky & Simon 2008). 
Bresolin (2011) has demonstrated that this scatter is not intrinsic but is instead due to the 
measuring errors of the very weak [OIII] A4363 lines needed to derive the electron tempera- 
ture in higher metallicity regions. Bresolin (2011) was able to measure this line in 8 regions 
inwards of 3 kpc, and derives the oxygen abundances shown in Figure [T^ by the dashed red 
line. He does show that if one ignores the direct electron temperatures, and instead uses 
the so-called -R23 method (Pagel et al. 1979) on a larger sample, one derives a much higher 
metallicity. However Bresolin (2011) emphasizes that the exact value depends upon which 
high-met allicity calibration of the R23 index one adopts, although he does find a nearly iden- 
tical slope as he does with the data for which there are direct temperature measurements. 
Magrini et al. (2010) argues that by restricting the analysis to HII regions with measurable 
[O HI] A4363 lines, one is automatically biasing the result towards lower metallicities, a 
selection effect which Bresolin (2011) acknowledges. Instead, Magrini et al. (2007) proposes 
adopting a two-component model (shown in Figure [T^ by the solid green lines), where the 
inner oxygen abundances are based upon the oxygen abundances of young stars. In Fig- 
ure [T3r right) we emphasize just the Magrini et al. (2007) and the Bresolin (2011) results, 
with vertical dashed lines to indicate the effective p values for the three sample regions. 
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We see from Figure [12] that the WC/WN ratio in the outer two bins (p = 0.38 and 
0.69) are similar to that seen in the LMC Indeed the metalhcities of these two bins are 
log(0/H) + 12 = 8.4 and 8.3, respectively, according to Figure [T3l perfectly consistent with 
the oxygen abundance of the LMC, 8.37 (Russell & Dopita 1990). However, the innermost 
bin (p = 0.16) must have a metallicity that is considerably greater than that of the LMC 
based on its high WC to WN ratio. Were we to linearly extrapolate from the SMC to the 
LMC we would expect a log(0/H) + 12=8.9-9.0. This is clearly at variance with that of the 
Bresolin (2011) result, but is nicely in accord with the two-component gradient proposed by 
Magrini et al. (2007). 

We can provide one other bit of information on the metallicity gradient of M33. It has 
long been known (Smith 1968c) that nearly all of the late-type WC stars are found inwards 
of the solar circle, while the early-type WCs are found outwards of the solar circle. Smith 
(1968c) supposes that this is due to the metallicity gradient in the Milky Way, as the result 
is consistent with the fact that the SMC and LMC (which are metal poor compared to the 
solar neighborhood) have only early-type WCs. Massey & Johnson (1998) note that it is only 
M31 (where the metallicity is higher) that contains an abundance of late-type WCs. The 
interpretation that the spectral subtype of the WC class is determined by initial metallicity 
has been given a theoretical underpinning by Crowther et al. (2002). Furthermore, late-type 
WCs have CIV A5806 lines that are weaker (smaller EWs) and skinnier (smaller FWHMs) 
than in early-type WCs, so a change in average spectral type can be demonstrated using a 
plot of these two quantities whether or not the spectral subtype has been well determined. 
Massey & Johnson (1998) used such a diagram to demonstrate that the WCs in the inner 
region of the M33 were later in type than those in the outer regions. With the considerably 
more complete data here, we revisit this issue in Figure [Til In the three upper panels we see 
a very clear progression from the inner regions (p < 0.25), where the WCs are dominated 
by stars with narrower CIV A5806 lines (i.e., mostly of later type), while the outer region 
(p > 0.50) are dominated by WCs with stronger and broader lines. The middle region 
(0.25 < p < 0.50) is intermediate, but more similar to the outer region than the inner. This 
provides a largely independent confirmation that the metallicity difference from the inner 
region to the outer region differs by some significant amount. It is hard to place an exact 
value on this, but we can compare the distributions to that of the Milky Way and LMC in 
the lower two panels. Clearly the distribution of the WCs in the inner region (upper left) 
is more like that of the Milky Way (lower left) than like that of the LMC (lower right), 
while the outer region of M33 (upper right) is more like that of the LMC. Thus, to the 
extent that the WC spectral properties are dominated by metallicity effects, the inner region 
of M33 (p < 0.25) must have an oxygen abundance that is more similar that of the solar 
neighborhood (log(0/H) + 12 = 8.7, Esteban & Peimbert 1995) than that of the LMC (8.4). 
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This provides additional, independent support for the higher metalhcity value for the center. 

In the following section we will therefore adopt the two-component oxygen gradient of 
Magrini et al. (2007). 

5.3. Comparison with the Evolutionary Models 

We can now compare our observed WC to WN ratios with the predictions of the Geneva 
evolutionary models. In this paper we have chosen to use only the galaxies where the WR 
content is well-known (i.e., SMC, LMC and now M33) when comparing the ratios with 
the models. For other galaxies in the Local Group, the surveys are thought be selectively 
incomplete for WNs (Massey & Holmes 2002). This is especially true for the Milky Way 
where the interstellar extinction further complicates the matter (Shara et al. 2009). Figure [T5] 
shows both the expected ratio of WCs to WNs as computed by Meynet & Maeder (2005) as 
well as what we observe, with the data taken from Table M 

As described by Meynet & Maeder (2005), the models fail for high WC to WN ratios 
(such as those believe to exist in M31 and the Milky Way), but there has always been the 
lingering doubt that incompleteness in the observational bias towards WCs might be possibly 
responsible (see discussion in, for example, Massey Sz Holmes 2003 and Meynet & Maeder 
2005). The result that the discrepancy exists in the inner part of M33 is new here. Note 
too that this result is robust to any uncertainty in the metallicity of the center of M33 (i.e.. 
Section [5^ as the rotating models fail to predict this high a value of the WC/WN ratio at 
any metallicity. As long as the WN content of the inner bin is relatively well known, as we 
argue above, then the problem exists. 

6. Summary and Conclusions 

M33 contains 206 confirmed WR stars, 55 of which were discovered as part of this survey. 
Based on our remaining candidates, we expect that there may be another half-dozen WRs 
that have yet to be found. This number is in accord with the ~190 WR stars we expect 
by scaling from the LMC based on luminosity. This suggests that the current massive star 
formation rates in M33 are comparable with those in the LMC per unit luminosity. Overall, 
the dominant WC subclass is WC4 while there are a few later type WCs (WC7s) located in 
the central region of M33. The WNs are dominated by early- type WNs, with some late- type 
WNs found predominantly in rich regions of massive stars, although a few are located in the 
field. 
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The WC/WN ratio changes dramatically with galactocentric distance in M33. This can 
be seen quite obviously in Figure [11] and in Table [51 We believe this is due to a metallicity 
gradient, consistent with the two-component model observed by Magrini et al. (2007). This 
is in accord with the basic premise of massive star evolutionary theory (the "Conti Scenario" , 
Conti 1976): in a high metallicity environment, mass-loss rates will be higher for a given 
luminosity and mass, and thus more stars will be found in the WC stage at a given time. 
Massey & Conti (1983) suggest the change in the WC/WN ratio with galactocentric distance 
might be due to a change in the IMF slope, but current studies in regions of star formation all 
suggest that the IMF slope is invariant at the upper end (Bastian et al. 2010). Additionally, 
WR stars can also be formed as a result of Roche-lobe mass transfer in a binary system, but 
we do not expect binary WRs to significantly influence M33's overall WC/WN ratio. 

While the Geneva models (Meynet & Maeder 2005) do an good job at predicting the 
WC/WN ratio at lower metallicities, they clearly fail at solar metallicities (i.e., log(0/H) + 
12 ~ 8.7) by not predicting a sufficiently high WC/WN ratio. Eldridge & Vink (2006) 
suggest that the problem with the evolutionary models is solved by adding a metallicity- 
dependence mass-loss rate to the stellar winds during the WR stage. However, Meynet et 
al. (2008) argue that such models severely underestimate the number of WRs relative to O- 
type stars, and instead proposes that the problem is alleviated if an LBV phase is included 
for the very massive stars that enter a WN phase while still H-burning. Including an LBV 
stage would then significantly shorten the predicted WN lifetime, increasing the predicted 
WC/WN ratio at higher metallicities. Meynet et al. (2008) report that tests of the Geneva 
models at solar metallicity confirm that this helps. 

In the future, we plan on conducting similarly complete surveys of M31 (log(0/H) + 12 ~ 
9.0), ICIO (log(0/H) + 12 = 8.25) and NGC6822 (log(0/H) + 12 = 8.25). This will allow 
us to further understand the relationship between metallicity and the WC to WN ratio and 
it will allow the theorists to improve heavily-relied-upon evolutionary models. 
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Section [^?T] and Susan Tokarz for reducing the spectra as described in Section [31 Additionally, 
we thank Carl Akerlof and Fang Yuan for their helpful suggestions on compiling and running 
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useful comments on this manuscript by Georges Meynet, Laura Magrini, Carol Kehrig, and 
Peter S. Conti. Howard Bond, Steve Howell, Bill Keel, and Scott Kenyon all made useful 
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CT 



WN Result 



Fig. 1. — Demonstration of the image subtraction technique. The first image shows a region 
of M33 through the CT filter, the second through the WN filter and the third is the result 
of using the image subtraction program. Three confirmed WN stars are outlined in red 
dashed circles. While the WR star on the lower left is clearly more prominent through the 
WN filter than through the CT filter, the other two WR stars appear by eye to be the 
same brightness. However, the image subtraction program is able to detect the stars' small 
brightness differences through the two different filters. Notice that while there are many 
stars in the field, only three show up prominently in the result image. 
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4000 5000 6000 7000 8000 

Wavelength (A) 

Fig. 2. — The spectra of two newly discovered quasars. Upper: the star 
J0133445.02+304928.0 has a redshift of ^ = 1.99, and its strong CIV A1550 hne had been 
shifted to 4645 A, placing it near the center of our WC detection filter. Lower: the star 
J013322.09+301651.4 has a redshift of ^ = 2.84, and its Lya line had been shifted to 4672 
A, placing it in the WN filter bandpass. That star had previously been called a WN-type 
Wolf-Rayet star (Massey & Conti 1983), with the Lya featured mistaken for a blend of NIII 
A4634,42 and He II A4686 by Massey & Conti (1983). 
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Fig. 3. — Two non-WR emission-line stars. Both stars show strong hnes of H and He I, 
with weak [Nil] AA6548, 6584 nebular emission. Left: J013423.00+304536.5 additionally 
shows He H emission hnes, and at My — —4.5 may be a symbiotic star. The counter- 
argument is that there is no sign of a late- type companion visible in the far red. Right: 
J013351. 46+304057.0 shows a broad feature that may include a He H A4686 component but 
not other He II lines. Its spectrum shows [SII], originating in a low density region. At 
My — —7.2 it is too luminous be a symbiotic star. 
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Fig. 4. — The spectra of representative WN stars. On the left we show the 3750-7400 A 
region of the spectrum, while on the right we show an expansion of just the blue region. 
The latter contains the N III A4634, 42, N IV A4058, and NV A4603, 19 hues which form the 
primary basis for the classification. We display one example for each of the WN subtypes 
found in M33. 
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Fig. 4. — Continued 
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Fig. 4. — Continued 



J013245.84+302019.5 
WC4+neb 





MM 



4000 4500 5000 5500 6000 
Wavelength (A) 



"S 10 



JO 13507.24+304500.9 

¥C5 



lvvJw■~^/^-^^j|^JV 



4000 4500 5000 5500 
Wavelength (A) 



_J013338. 20+303112. e 

_WC6 




4000 4500 5000 
Wavelength 



5500 6000 



(A) 



Fig. 5. — Three WC stars. The classification of WC stars is based upon the relative 
strengths of C III A5696, O V A5592, and C IV A5806. Of these three, J013338.20+303112.8 
(WC6) was previously known (MC 35). The unmarked narrow lines in the spectrum of 
J013245. 84+302019. 5 are nebular. We display examples of each WC subtype found. 
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Fig. 6. — The spectrum of the 04 If star J013333. 58+304219. 2, previously (mis) classified as 
a WN star by Massey & Johnson (1998). Although superficially one might be inclined to 
describe the spectrum (left) as that of a late-type WN star, the presence of He I and He II 
absorption features (right) shows that this is actually an 0-type star, which we classify here 
as 04 If. This is the earliest type supergiant classified to date in M33. 
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Fig. 7. — Three Ofpe/WN9 stars. The Ofpe/WN9 class is marked by properties intermedi- 
ate between that of Of-type stars and very late-type WN type. J013509. 73+304157. 3 was 
previously classified by Massey et al. (2007b) as Ofpe/WN9; the object is sometimes referred 
to as "Romano's star" based on the discovery spectrum described by Romano (1978). The 
other two stars whose spectra are shown here are newly discovered here. 
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Fig. 8.— Changes in the spectrum of J013245.41+303858.3. At the top our new (2010) 
spectrum is shown, while at the bottom the old (1995) spectrum of Massey et al. (1996) is 
shown. The fact that Balmer hydrogen lines are weaker in the new spectrum than in the old 
may be just an artifact of sky subtraction being non-local with fiber instruments, and the 
2x larger diameter fibers used in the 1995 observation. But as the figure on the right shows, 
the He I A4471 line has also gone from being P Cygni to pure absorption. 
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Fig. 9. — The WN2pec star J013510. 27+304522. 9. This star is very unusual in having strong 
O VI A3811, usually an indication of very high temperature and oxygen abundance. 
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Fig. 10. — Sensitivity of Current Survey. The log of the flux in the He II A4686 emission 
line is plotted against absolute visual magnitude My- The line flux has been approximated 
by log(— EW) — Mv/2.5. The open circles denote the WNs detected in our survey (both 
new and previously known). The solid circles are the Of-type stars. The data for the 
SMC (red asterisks) and LMC (blue x's) are shown for comparison, where the data comes 
from Massey et al. (2003) and Conti & Massey (1989), respectively. The current survey has 
detected WNs with line fluxes as weak or weaker than WNs known in the Magellanic Clouds, 
arguing that the survey has sufficient sensitivity to allow a meaningful determination of the 
relative number of WC and WN stars. 



Fig. 11. — Location of known WC and WN stars in M33. WN stars are represented as blue 
xs while WC stars are represented as red +s. The green ovals represent distances of p = 0.25 
(1.9 kpc) and p = 0.50 (3.8 kpc) within the plane of M33. Our survey extended well outside 
the bounds of this figure, except for the small missing region in the lower left corner. 
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Fig. 12. — WC/WN ratio vs galactoccntric distance. The ratio of the number of WC-type 
and WN-type WRs is shown plotted against the distance from the center within the plane 
of M33. A value of p = 1 corresponds to a distance of 7.53 kpc. The error bars are based 
on a 5% incompleteness. The values of the WC/WN ratios for the LMC and SMC are also 
shown. Note that the LMC and SMC have log(0/H) + 12 = 8.37 and 8.13, respectively, and 
thus this plot implies that the metallicity of the two outmost bins must be similar to that 
of the LMC, while that of the inner-most bin is significantly (~ 0.5 dex) greater than that 
of the LMC. 
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Fig. 13. — Oxygen abundances vs galactocentric distance. Left: The oxygen abundance 
gradient from Crockett et al. (2006) is shown by a green dotted hne, that of Rosolowsky 
& Simon (2008) is shown by a short-dashed purple hne, that of Magrini et ah (2007) with 
a sohd blue hne, and that of Bresolin (2011) by a long-dashed red line. Right: We show 
the oxygen abundance gradient from Magrini et al. (2007) and Bresolin (2011), with the 
vertical dash lines indicating the average p values of the three regions denoted in Figure 
and Table |9l 
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Fig. 14. — WC line strength vs line width. The line strength of the CIV A5806 line is plotted 
against its line width with the spectral subtypes indicated. The upper panels show the data 
for the three M33 regions (either based upon measurements here or from the references 
indicated in Tabled]), while the while two panels show the data for the Milky Way and LMC 
(data from Conti & Massey 1989, tables 3 and 4). The distribution for the inner region of 
M33 (p < 0.25, upper left) is more similar to that of the Milky Way (lower left) than the 
LMC (lower right), while the outer region (upper right) are more like that of the LMC. 
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Fig. 15. — The WC/WN ratio versus metallicity. The points are from Tabled and the hne 
shows the evolutionary model results from Meynet & Maeder (2005) for an initial rotation 
of vsini = 300 km s~^. 



Table 1. Observed Stars and Success Rates 



All B-V <1 



Priority^ 


Total # 


Observed 


Total # 


Observed 


Winners'^ 


Losers 


% Winners 


# WRs 


% WRs 


1 


92 


83 


83 


77 


68 


9 


88.3 


55 


68.8 


2 


104 


66 


93 


56 


6 


49 


10.7 


2 


3.6 


3 


105 


64 


82 


45 


1 


44 


2.2 





0.0 


4 


75 


25 


67 


19 


2 


17 


10.5 





0.0 


5 


132 


30 


132 


30 


30 





100.0 


29 


96.7 



^1-found both by image subtraction and by photometry. 2-High significance (> 50") photometry candidate. 3 — 
Found by image subtraction and not by photometry. 4 — Lower significance (5-6cr) photometry candidate. 5 — Known 
WRs. 

•"'Winner" means only that there is emission within the emission- fine bandpass, not that the star is a WR. 
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Table 2. Non-Of and Non-WR stars with Emission 



Star 


V 


B-V 


U-B 


Type 


Comment 


J013322.09+301651.4 


19.42 


0.32 


0.03 


QSO 


z=2.84 MC25 


J013331.23+303334.2 


20.25 


0.40 


-1.87 


HIT region 


Nebular He II A4686 


J013351.46+304057.0 


17.73 


0.08 


-0.93 


Unknown 


See text 


J013423.00+304536.5 


20.47 


0.27 


-1.06 


Symbiotic? 


See text 


J013445.02+304928.0 


20.11 


0.21 


-1.04 


QSO 


z=1.99 


J013459.47+303701.9 


18.37 


0.21 


-1.23 


Hot LBV cand 


Massey et al. (2007b) 



Table 3. New WR and Of-type Spectroscopic Results 



Star 


Type 


Previous 




cm and/or He II 


CIV A5806 


Comments 






Type 


ID 


Ref^ 


log(-EW) 


FWHM 


log(-EW) 


FWHM 




J013232.07+303522.4 


WN2 








1.7 


37 








J013232.13+303514.3 


WN3+neb 








1.8 


36 


1.5 


50 




J013233.24+302652.3 


WN2+abs 








1.0 


30 








J013241. 39+303416.2 


WN3h 








2.1 


26 


0.4 


25 


H-rich 


J013245.41+303858.3 


Ofpe 


Ofpe/WN9 


UIT 008 




0.3 


14 








J013245.74+303854.4 


WN2+07 


WN 


MC 3 




1.1 


35 








J013245.84+302019.5 


WC4+neb - 








3.0 


76 


3.5 


65 


OV present 


J013251.04+303535.1 


Olf 








0.5 


7 






Neb subtraction poor 


J013255.91+303155.9 


WC4 


wc 


MC 4 




2.3 


70 


2.5 


110 


OV present ^ 

1 


J013256.35+303535.4 


WN6+abs 


WN6+abs 


OB 21-65 




1.1 


14 


0.0 


19 


J013257.32+304418.1 


WN3+neb 








2.1 


45 






1 


J013257.88+303549.8 


WC4 


WC4-5 


MC 6 




1.9 


72 


2.5 


100 


OV present 


J013302.28+301119.1 


WN3+neb 


WN3-4 


KlO-01 


1 










Mai w Hell neb emis 


J013302.67+303120.1 


WN3 


WNE 


MC 7 




2.3 


37 


1.0 


34 




J013303.21+303408.6 


WN6 


WNL 


MC 10 




1.7 


15 


0.8 


19 




J013307.50+304258.5 


WN2 


WNE 


UIT 041 


2 


1.1 


62 






Neb. Balmer but little [OIII] 


J013307.68+303315.4 


WN4 


WN 


MC 14 




1.9 


26 


1.0 


32 




J013310.77+302734.2 


WN4 








1.8 


23 


1.0 


44 




J013311.13+303916.8 


WN3+abs 








0.6 


22 








J013311.29+303146.9 


WN4+B2I 








0.8 


23 








J013312.44+303848.0 


WN3+neb 


WN4 


N592-WR2 


3 


1.5 


35 


0.5 


43 


NIV present 


J013312.54+303900.3 


WN/C+abs 








1.3 


24 


1.2 


37 




J013312.61+304531.0 


WN4 








1.7 


24 


1.4 


34 




J013314.56+305319.6 


WN3+neb 








1.5 


23 









Table 3 — Continued 



Star 


Type 


Previous 




cm and/or He II 


CIV A5806 




Comments 






Type 


ID 


ReP 


log(-EW) 


FWHM 


log(-EW) FWHM 




J013315.02+303906.9 


WC5-6 








Z. 1 


DO 


2.4 


62 




J013315.43+302300.1 


WN3 








1 
l.Z 










J013315.55+304514.1 


WN3+ncb 


W IN 






z.u 


OQ 
Z6 








J013315.82+305644.8 


WN4.5+0 


vviN4.o-t-neD 






1 
l.Z 


/in 


0.1 


30 


Strong neb 


J013316.17+304751.8 


WC4 




A/TP 09 




Q 


fiO 
OZ 


3.1 


67 




J013320.53+303201.5 


earlyOf 
















NIII strong 


J013324.04+305030.7 


WN8 








1 
l.Z 


1 n 








J013326.60+303550.3 


WN6 








1 c; 

i.O 


1 
IZ 


0.4 


19 




J013329.88+303147.3 


09:If+neb 








n /I 


1 

iO 








J013330.36+303128.6 


WN4 








1 7 


O/l 
Z4 








J013332. 23+304131. 2 


061(f) 
















NIII strong 


J013333.58+304219.2 


04If 


WIN 


iViJ uzz 




U.o 


7 
( 








J013334.77+304101.5 


09If 
















CIV 0.5EW abs 


J013334.93+310042.8 


WN4 








1 c; 

1.0 


00 

zz 


0.4 


29 




J013335.23+310037.6 


WN3 








2.3 


27 


0.8 


27 




J013336.67+304302.6 


WN4+0 








1.1 


24 


0.3 


40 




J013337.34+303527.1 


WN4 


WN4 


MCA 5 


4 


1.7 


25 


0.6 


42 




J013337.81+302831.9 


WN4 








1.8 


25 


1.0 


35 




J013338.20+303112.8 


WC6 


WC6 


MC 35 


4 


2.4 


37 


2.1 


36 




J013340.19+303134.5 


WC4 


WC4 


MC 42 


4 


2.4 


84 


2.3 


70 


OV present 


J013340.32+304600.9 


WN3+neb 


WN? 


MJ B13 




2.0 


47 








J013343.33+304450.5 


WN4 


WN4 


MCA 8 


4 


1.7 


25 


0.2 


19 




J013350.71+305636.7 


WN3+neb 








2.4 


55 


1.6 


93 




J013352.43+304351.7 


WN8 


WNL 


MCA 12 




1.0 


11 









Table 3 — Continued 



Star 


Type 


Previous 


cm and/or He II 


CIV A5806 




Comments 






Type ID ReP 


log(-EW) 


FWHM 


log(-EW) FWHM 




J013355.33+302001.0 


WN3+neb 




A 




1.4 


48 




J013355.87+304528.4 


WN7 


WNL+neb ■ ■ ■ 2 


1 n 
i.U 


1 Q 
lo 








J013355.94+302732.1 


WN3+B2I 




1.1 


ZD 


0.3 


37 




J013356.33+303420.8 


WN3 




1.1 










J013359.60+303435.4 


WN5 




1.6 


oo 
zz 


0.7 


27 




J013359.79+305150.0 


WN2+09I 




U.o 


oU 








J013401.92+303819.3 


05:If+neb 




U.o 


y 








J013402.93+305126.2 


WN3+neb 




1 
Z.l 


1 


1.3 


84 




J013406.72+304154.5 


06If 




-U.l 


Q 






CIVO.SEWabs I 


J013406.80+304727.0 


WN7+neb 


WN7+neb UIT 303 • ■ • 


1.1 


01 
Zl 






Could be Of 


J013408.23+305234.2 


WN4 




1 

z.l 


OU 


1.2 


39 


1 
1 


J013408.90+304732.0 


WN3 




z.u 


/I ^ 


1.2 


57 




J013410.72+305240.7 


WN3 




1 

z.l 


OO 


1.4 


47 




J01341 1.14+304637.5 


WC4 




1 

z.l 


00 


2.1 


45 


No CHI of OV 


J013413.06+305230.0 


05If 












Hell A4686 lA EW in abs 


J013415.85+305522.9 


WN4 




1.7 


23 


0.7 


28 




J013416.07+303642.1 


Ofpc/WN9 














J013416.28+303646.4 


WC6+neb 


WC5 MC 70 4 


2.6 


45 


2.3 


34 




J013417.21+303334.7 


WN3+neb 




2.0 


29 


1.2 


42 


NIV present 


J013418.37+303837.0 


Ofpe/WN9 














J013418.44+303430.6 


05If 




0.5 


8 








J013419.16+303127.7 


WN3 




2.0 


27 








J013419.58+303801.5 


WN7 




0.8 


8 






Not Of as other Hell emis 


J013419.68+303343.0 


WN3-4+neb 




1.8 


43 






NIV present 



Table 3 — Continued 



Star 


Type 




Previous 


cm and/or He II 


CIV A5806 




Comments 






Type 


ID ReP 


log(-EW) 


FWHM 


log(-EW) FWHM 




J013421.21+303758.2 


WN5+neb 






z.u 


ZO 


1.3 


32 




J013421.97+303314.3 


WN4 






i.O 


0/1 


1.0 


33 




J013422.37+303313.7 


WN3+abs 


WN4 


KlO-03 1 


1 

i.z 




0.3 


38 


Hel absorption 


J013423.02+304650.0 


WN4 






i.O 


ZD 


0.3 


40 




J013425.11+301950.3 


WN3 






1 Q 

i.y 


Zo 






No NIV 


J013426.96+305256.7 


WN4.5 






1 A 


01 

zi 


0.1 


15 




J013427.08+305315.7 


08If 






-U.4 











J013427.30+305229.2 


WN3+neb 






1 


Ofi 
ZD 


1.3 


37 


NIV present 


J013429.56+304145.0 


WC4 






z.o 


DD 


2.9 


73 


CHI present 1 


J013431.45+305716.7 


WN3+neb 


WN 


MC 73 


Q 
Z.O 


40 


1.5 


74 


J013432.18+304903.7 


WC4 






A 


oO 


3.0 


140 


1 
1 


J013432.68+304655.4 


Olf+neb 


WN6 


N604-WR10 3 


U. / 


1 A 
iO 


0.4 


23 




J013433.22+310019.3 


WN7+neb 






i.O 


1 1 
ii 


0.0 


18 


H-rich 


J013433.82+304656.3 


WC4 


WC4 


N604-WR7 3 


z.u 


7Q 
/ O 


2.0 


50 


No CHI or OV 


J013435.15+304705.1 


08If 






-0.3 


8 






NIII strong 


J013438.18+304953.2 


WN5 






2.1 


23 


1.4 


34 




J013438.98+304119.8 


WC4 






2.1 


48 


2.6 


95 


OV present 


J013439.32+304408.6 


07If 






-0.4 


8 






NIII strong 


J013439.97+304749.8 


04111(f) 














Hell A4686 abs but NIII emis 


J013440.42+304321.9 


WN4+neb 


WN 


MCA 16 


1.6 


22 


0.9 


32 




J013443.51+304919.4 


WN4 






2.0 


24 








J013444.28+303757.2 


WN3 






1.6 


25 








J013444.61+304445.4 


WC5-6 


wc 


MC 77 


2.9 


55 


2.8 


43 




J013447.32+310748.8 


WN3 






1.8 


25 









Table 3 — Continued 



Star 


Type 




Previous 


cm and/or He II 


CIV A5806 


Comments 


Type 


ID Ref^ 


log(-EW) 


FWHM 


log(-EW) FWHM 


J013448.68+302819.8 


09If 






0.0 


24 






J013507.24+304500.9 


WC5 






3.1 


55 


2.9 43 


CHI and OV present 


J013509.73+304157.3 


Ofpe/WN9 


Ofpe/WN9 


Romano's Star 1 










J013510.27+304522.9 


WN2pec 






2.3 


50 




Very strong VI A3811 



'^References for old spectral types are from Massey & Johnson 1998 unless otherwise noted. 1 — Kehrig et al. 2011; 2 — Massey et al. 2007b; 
3— Drissen et al. 2008; 4— Abbott et al. 2004. 

I 
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Table 4. Stars No Longer Considered WRs 



Star 


Old 




New 




Other ID 




Type 


Ref.^ 


Type 


Ref.^ 




J013322.09+301651.4 


WN 


1 


QSO 


2 


[MJ98] 37, MC25 


J013333.58+304219.2 


WN 


3 


04 If 


2 


[MJ98] 46, C22 


J013344.68+304436.7 


WN8 


4 


If 


5 


[MJ98] 77, OB66-25 


J013416.07+303642.1 


WNll 


6 


P Cyg LBV 


7 


B517, H 108 


J013432.68+304655.4 


WN6 


8 


Olf 


2 


N604-WR10 



^1— Massey & Conti 1983; 2— This work; 3— Massey & Johnson 1998; 4— Massey et al. 
1996; Abbott et al. 2004; 6— Crowther et al. 1997; 7— Massey et al. 2007b; 8— Drissen et al. 
2008. 



Table 5. WR Stars in M33 



Star 




Type 


Ref'^ 










OB^ 


MJg 


ID*^ Comment 


J013232.07+303522.4 


0.88 


WN2 


1 


20.99 


20.76 


-4.0 


-4.3 


131 






J013232.13+303514.3 


0.88 


WN3+neb 


1 


20.74 


20.43 


-4.2 


-4.6 


131 






J013233.24+302652.3 


0.85 


WN2+abs 


1 


21.22 


22.23 


-3.8 


-2.8 


123 






J013237.62+304004.3 


0.87 


WC 


2 


19.77 


20.00 


-5.2 


-5.0 


Fid 


WRl 


MCA 1 


J013237.72+304005.7 


0.87 


Ofpe/WN9 


2 


17.75 


20.00 


-7.2 


-5.0 


132 


WR2 


UIT 003 


J013240.82+302454.3 


0.79 


WNE 


2 


20.72 


20.51 


-4.3 


-4.5 


121 


WR3 


MC 1 


J013241. 39+303416.2 


0.76 


WN3h 


1 


20.75 


20.41 


-4.2 


-4.6 


Fid 






J013241.95+304024.6 


0.83 


WCE 


2 


20.55 


21.09 


-4.4 


-4.0 


132 


WR4 


MC 2 


J013245.41+303858.3 


0.76 


Ofpe/WN9 


2 


17.61 


17.29 


-7.4 


-7.8 


27 


WR5 


UIT 008 


J013245. 74+303854.4 


0.76 


WN2+07 


2 


18.00 


17.74 


-7.0 


-7.3 


27 


WR6 


MC 3 


J013245.84+302019.5 


0.82 


WC4+neb 


1 


20.32 


21.11 


-4.7 


-3.9 


(118) 






J013255.91+303155.9 


0.60 


WC4 


1 


21.69 


21.61 


-3.3 


-3.4 


129 


WR7 


MC 4 


J013256.35+303535.4 


0.60 


WN6+abs 


1 


18.93 


18.90 


-6.0 


-6.1 


21 


WR8 


OB 21-65 


J013256.84+302724.9 


0.62 


WN 


2 


21.34 


21.07 


-3.6 


-4.0 


122 


WR9 


MC 5 


J013257.32+304418.1 


0.72 


WN3+neb 


1 


21.07 


20.97 


-3.9 


-4.1 


Fid 






J013257.88+303549.8 


0.58 


WC4 


1 


19.55 


19.44 


-5.4 


-5.6 


21 


WRIO 


MC 6 


J013300.20+303015.3 


0.56 


WN8 


2 


18.78 


18.49 


-6.2 


-6.6 


128 


WRll 


MJ El 


J013302. 28+301119.1 


0.98 


WN3+neb 


1 


19.87 


19.50 


-5.1 


-5.5 


MAI 




KlO-01 


J013302.67+303120.1 


0.53 


WN3 


2 


20.79 


20.41 


-4.2 


-4.6 


128 


WRl 2 


MC 7 


J013302.94+301122.8 


0.98 


WNE 


2 


16.14 


20.41 


-8.8 


-4.6 


MAI 


WR14 


MC8 


J013303.19+301124.2 


0.98 


WNE 


2 


17.23 


20.41 


-7.7 


-4.6 


MAI 


WR15 


MC9 


J013303.21+303408.6 


0.52 


WN6 


2 


20.89 


20.68 


-4.1 


-4.4 


20 


WR13 


MC 10 


J013303.71+302326.3 


0.64 


WC 


2 


21.75 


21.45 


-3.2 


-3.6 


Fid 


WR16 


MC 11 


J013304.98+303159.8 


0.50 


WNE 


2 


21.79 


21.45 


-3.2 


-3.6 


19 


WRl 7 


MC 12 


J013305.67+302857.6 


0.52 


WNE 


2 


20.41 


20.30 


-4.6 


-4.7 


Fid 


WR18 


MC 13 



Table 5 — Continued 



Star 




Type 


Ref^ 








Ma4750^ 


OB^ 




ID'' 


Comment 


J013307.50+304258.5 


0.58 


WN2 


2 


17.33 


20.11 


-7.6 


-4.9 


28 


WR19 


UIT 041 




J013307.68+303315.4 


0.47 


WN4 


2 


21.35 


21.06 


-3.6 


-4.0 


19 


WR20 


MC 14 




J013307.80+302951.1 


0.49 


WNE 


2 


20.88 


20.60 


-4.1 


-4.4 


127 


WR21 


MC 15 




J013308.56+302805.6 


0.51 


WNE 


2 


20.68 


20.40 


-4.3 


-4.6 


Fid 


WR23 


MC 16 




J013309.14+304954.5 


0.74 


Ofpe/WN9 


2 


17.91 


20.40 


-7.1 


-4.6 


28 


WR22 


UIT 045 




J013310.74+303900.2 


0.47 


WN 


2 


20.89 


20.63 


-4.1 


-4.4 


Fid 


WR24 


MCA 2 




J013310.77+302734.2 


0.51 


WN4 


1 


20.98 


20.82 


-4.0 


-4.2 


115 


WR144 


MJ-E16 


Prev. WN cand. 


J013311.13+303916.8 


0.46 


WN3+abs 


1 


21.17 


20.82 


-3.8 


-4.2 


59 








J013311.29+303146.9 


0.44 


WN4+B2I 


1 


18.73 


18.39 


-6.2 


-6.7 


17 


WR145 


MJ-E34 


Prev. WN cand. 


J013311.44+304856.9 


0.69 


WNE+B3 I 


2 


16.49 


18.39 


-8.5 


-6.7 


135 


WR26 


UIT 052 




J013311.85+303852.7 


0.45 


WNL 


2 


18.14 


18.00 


-6.8 


-7.0 


59 


WR25 


MC 17 


N592-WR1 


J013312. 21+302740.5 


0.49 


WN 


2 


21.14 


20.90 


-3.8 


-4.1 


115 


WR27 


MJ E5 




J013312.44+303848.0 


0.44 


WN3+neb 


5 


19.22 


18.76 


-5.8 


-6.3 


59 




N592-WR2 




J013312.54+303900.3 


0.44 


WN/C+abs 


1 


19.78 


18.76 


-5.2 


-6.3 


59 








J013312.61+304531.0 


0.58 


WN4 


1 


20.44 


20.10 


-4.5 


-4.9 


29 








J013312.95+304459.4 


0.56 


WN 


2 


21.31 


21.08 


-3.7 


-4.0 


29 


WR28 


MC 18 




J013314.34+302955.3 


0.44 


WN 


2 


20.22 


20.04 


-4.8 


-5.0 


127 


WR29 


MJ E6 




J013314.56+305319.6 


0.79 


WN3+neb 


1 


19.92 


19.88 


-5.1 


-5.2 


137 








J013315.02+303906.9 


0.42 


WC5-6 


1 


19.78 


19.96 


-5.2 


-5.1 


59 


WR30 


MC 19 




J013315.31+304503.5 


0.54 


WN3, WN/CE 


2 


18.53 


18.23 


-6.4 


-6.8 


29 


WR31 


MC 20 




J013315.43+302300.1 


0.59 


WN3 


1 


22.08 


18.23 


-2.9 


-6.8 


116 








J013315.55+304514.1 


0.54 


WN3+neb 


1 


21.18 


21.06 


-3.8 


-4.0 


29 


WR32 


MC 21 




J013315.82+305644.8 


0.89 


WN4.5+0 


1 


16.78 


16.60 


-8.2 


-8.4 


IC 


WR33 


MC 23 




J013316.17+304751.8 


0.61 


WC4 


1 


20.46 


21.13 


-4.5 


-3.9 


30 


WR34 


MC 22 




J013316.48+303221.3 


0.39 


WNL 


2 


20.19 


20.09 


-4.8 


-5.0 


17 


WR35 


MCA 3 





Table 5 — Continued 



Star 




Type 


Ref^ 














ID^ 


Comment 


J013318.50+302658.2 


0.47 


WC4-5 


2 


19.49 


19.31 


-5.5 


-5.7 


114 


WR36 


MC 24 




J013324.04+305030.7 


0.61 


WN8 


1 


19.07 


18.94 


-5.9 


-6.1 


136 








J013326.60+303550.3 


0.26 


WN6 


1 


20.11 


19.96 


-4.9 


-5.1 


(54) 








J013326.67+304040.3 


0.30 


WC6 


6 


21.14 


21.50 


-3.8 


-3.5 


61 


WR38 


AM 1 




J013327.26+303909.1 


0.27 


Ofpc/WN9 


2 


17.95 


17.80 


-7.0 


-7.2 


58 


WR39 


MJ C7 




J013327.76+303150.9 


0.31 


WN 


2 


19.72 


19.15 


-5.3 


-5.9 


15 


WR40 


MC 27 




J013330.36+303128.6 


0.30 


WN4 


1 


20.94 


20.62 


-4.0 


-4.4 


12 








J013332.64+304127.2 


0.25 


WNL/LBV 


2 


18.99 


18.89 


-6.0 


-6.2 


62 


WR41 


AM 2 


N595-WR6 


J013332.82+304146.0 


0.26 


WN7 


5 


18.20 


18.01 


-6.8 


-7.0 


62 


WR42 


AM 3 


N595-WR5 


J013332.97+304136.1 


0.25 


WNL 


2 


18.15 


18.00 


-6.8 


-7.0 


62 


WR43 


AM 4 


N595-WR4 


J013333.31+304129.7 


0.24 


WC6 


5 


19.34 


19.64 


-5.6 


-5.4 


62 


WR45 


AM 5 


N595-WR3 


J013333.49+304133.4 


0.24 


WNL? 


2 


19.81 


19.64 


-5.2 


-5.4 


62 


WR44 


MC 30 




J013333.75+304134.0 


0.24 


WNLx2? 


5 


17.42 


19.64 


-7.6 


-5.4 


62 


WR47 


MC 31 


N595-WR2ab 


J013333.81+304129.6 


0.24 


WN7 


5 


19.03 


19.64 


-5.9 


-5.4 


62 




N595-WR9 




J013334.04+304117.2 


0.23 


WN 


2 


19.78 


19.87 


-5.2 


-5.2 


62 


WR48 


MCA 4 


N595-WR8 


J013334.26+304138.1 


0.24 


WN7 


5 


18.38 


19.87 


-6.6 


-5.2 


62 


WR49 


AM 6 


N595-WR1 


J013334.27+304136.7 


0.24 


WNL+neb 


4 


16.97 


19.87 


-8.0 


-5.2 


62 : 








J013334.28+303347.5 


0.23 


WN5 


6 


21.08 


21.22 


-3.9 


-3.8 


14 


WR50 


MJ G13 




J013334.31+304130.4 


0.23 


WN8 


3 


19.52 


19.84 


-5.5 


-5.2 


62 


WR51 


AM 7 


N595-WR7 


J013334.93+310042.8 


0.88 


WN4 


1 


19.04 


19.32 


-5.9 


-5.7 


140 








J013335.23+310037.6 


0.88 


WN3 


1 


20.88 


21.06 


-4.1 


-4.0 


140 








J013335.47+304220.1 


0.24 


WC4 


6 


20.23 


20.31 


-4.7 


-4.7 


62 


WR52 


AM 8 




J013335.73+303629.1 


0.17 


wc 


2 


18.83 


20.31 


-6.1 


-4.7 


Fid 


WR53 


MC 34 




J013336.67+304302.6 


0.25 


WN4+0 


1 


19.27 


19.13 


-5.7 


-5.9 


(63) 








J013337.34+303527.1 


0.17 


WN4 


1 


21.11 


21.03 


-3.9 


-4.0 


50 


WR54 


MCA 5 





Table 5 — Continued 



Star 




Type 


Ref^ 


V 






Ma4750^ 






ID** Comment 


J013337.81+302831.9 


0.38 


WN4 


1 


21.34 


21.12 


-3.6 


-3.9 


Fid 






J013338.20+303112.8 


0.29 


WC6 


1 


20.65 


20.90 


-4.3 


-4.1 


12 


WR55 


MC 35 


J013339.30+303554.9 


0.15 


WN4 


6 


20.08 


19.87 


-4.9 


-5.2 


50 


WR56 


MC 36 


J013339.52+304540.5 


0.31 


B0.5Ia+WNE 


2 


17.50 


19.87 


-7.5 


-5.2 


Fid 


WR57 


UIT154 


J013339.70+302101.2 


0.66 


WNE 


2 


20.34 


20.30 


-4.6 


-4.7 


112 


WR58 


MCA 6 


J013339.95+303138.5 


0.27 


WN4b 


6 


19.67 


19.11 


-5.3 


-5.9 


12 


WR59 


MC 37 


J013340.04+303121.3 


0.28 


WNE+abs 


6 


19.76 


19.54 


-5.2 


-5.5 


12 


WR60 


MJ G3 


J013340.07+304238.4 


0.20 


WC5 


6 


20.00 


20.01 


-5.0 


-5.0 


Fid 


WR61 


AM 9 


J013340.19+303134.5 


0.28 


WC4 


6 


19.90 


19.87 


-5.1 


-5.2 


12 


WR62 


MC 42 


J013340.21+303551.8 


0.15 


WC4-5 


6 


19.34 


19.12 


-5.6 


-5.9 


50 


WR63 


MC 40 


J013340. 23+304102.0 


0.16 


WN 


2 


20.74 


20.66 


-4.2 


-4.4 


Fid 


WR64 


MC 38 


J013340.28+304053.4 


0.15 


wc 


2 


20.76 


21.34 


-4.2 


-3.7 


Fid 


WR65 


AM 10 


J013340.32+304600.9 


0.31 


WN3+neb 


2 


21.74 


22.06 


-3.2 


-3.0 


Fid 


WR66 


MJ B13 


J013340.69+304253.7 


0.20 


WN 


2 


20.56 


20.58 


-4.4 


-4.5 


Fid 


WR68 


AM 11 


J013341.65+303855.2 


0.10 


WN 


2 


20.16 


20.08 


-4.8 


-5.0 


Fid 


WR67 


MJ C3 


J013341.83+304154.8 


0.16 


WC6-7 


2 


20.56 


20.95 


-4.4 


-4.1 


Fid 


WR69 


AM 13 


J013341.91+304202.7 


0.17 


WN 


2 


20.30 


20.50 


-4.7 


-4.5 


Fid 


WR70 


MCA 7 


J013342.53+303314.7 


0.22 


WC4-5 


2 


18.99 


18.72 


-6.0 


-6.3 


9 


WR71 


MC 44 


J013343.19+303906.2 


0.09 


WN 7+abs 


2 


16.85 


18.30 


-8.1 


-6.7 


Fid 


WR72 


W91 129 


J013343.21+303900.5 


0.08 


WC5 


6 


18.63 


18.58 


-6.3 


-6.5 


Fid 


WR73 


MJ C4 


J013343.33+304450.5 


0.25 


WN4 


6 


21.08 


21.02 


-3.9 


-4.0 


66b 


WR74 


MCA 8 


J013343.34+303534.1 


0.14 


WN4.5+06-9 


4 


17.48 


21.02 


-7.5 


-4.0 


50 


WR75 


UIT 177 


J013344.40+303845.9 


0.07 


WC4 


6 


20.40 


20.36 


-4.6 


-4.7 


Fid 


WR76 


MC 45 


J013345.24+303841.1 


0.06 


WN7 


6 


20.60 


20.57 


-4.4 


-4.5 


Fid 


WR78 


MCA 9 


J013345.58+303451.9 


0.16 


WNL 


2 


20.75 


20.54 


-4.2 


-4.5 


49 


WR79 


MCA 10 



Table 5 — Continued 



Star 




Type 


Ref^ 








Ma4750^ 


OB^ 




ID** Comment 


J013345.99+303602.7 


0.12 


WN4b 


6 


20.43 


20.14 


-4.5 


-4.9 


48 


WR80 


MC 46 


J013346.20+303436.5 


0.17 


WNE 


2 


21.24 


21.28 


-3.7 


-3.8 


49 


WR81 


MCA 11 


J013346.55+303700.5 


0.09 


WC6 


6 


20.81 


20.94 


-4.2 


-4.1 


48 


WR82 


MC 47 


J013346.80+303334.5 


0.21 


WN6/C4 


6 


20.06 


20.09 


-4.9 


-5.0 


8 


WR83 


MC 48 


J013347.15+303702.8 


0.09 


WC6 


6 


20.75 


20.89 


-4.2 


-4.2 


48 


WR84 


MC 49 


J013347.67+304351.3 


0.18 


WN6 


2 


18.55 


18.41 


-6.4 


-6.6 


66 


WR85 


OB 66 F-61 


J013347.83+303338.1 


0.21 


WN4 


6 


20.01 


19.70 


-5.0 


-5.3 


8 


WR86 


MJ 08 


J013347.96+304506.6 


0.22 


WN 


2 


21.26 


21.38 


-3.7 


-3.7 


66c 


WR87 


MC 51 


J013348.85+303949.6 


0.03 


WN 


2 


19.62 


21.38 


-5.4 


-3.7 


Fid 


WR88 


MJ X4 


J013350.08+303818.9 


0.05 


WN+08-9 


2 


18.55 


18.33 


-6.4 


-6.7 


142 


WR90 


MJ X6 


J013350.08+303856.2 


0.02 


WC7+abs 


6 


19.91 


19.95 


-5.1 


-5.1 


142 


WR89 


MC 52 


J013350.23+303342.4 


0.22 


WN8-9 


6 


21.14 


21.92 


-3.8 


-3.1 


Fid 


WR91 


MJ 09 


J013350.26+304134.7 


0.08 


WC5-6 


6 


19.45 


19.36 


-5.5 


-5.7 


64 


WR92 


MC 53 


J013350.42+303834.2 


0.04 


Bl la+WN 


6 


17.95 


16.60 


-7.0 


-8.4 


142 


WR93 


UIT 213 


J013350.71+305636.7 


0.64 


WN3+neb 


1 


20.05 


19.63 


-4.9 


-5.4 


138 






J013351.28+303811.8 


0.06 


WC7+abs 


6 


19.77 


19.03 


-5.2 


-6.0 


142 


WR94 


MC 54 


J013351.84+303328.4 


0.24 


WC6 


6 


20.88 


21.03 


-4.1 


-4.0 


Fid 


WR95 


MC 55 


J013352.01+304023.5 


0.03 


WC4+abs 


6 


21.39 


21.93 


-3.6 


-3.1 


64 


WR96 


MC 56 


J013352.43+304351.7 


0.15 


WN8 


6 


19.08 


19.08 


-5.9 


-6.0 


65 


WR97 


MCA 12 


J013352.71+303907.3 


0.03 


WNE 


2 


20.72 


21.05 


-4.3 


-4.0 


142 


WR98 


MC 50 


J013352.71+304502.0 


0.20 


WC4 


6 


20.53 


20.70 


-4.4 


-4.3 


Fid 


WR99 


MC 57 


J013352. 75+304444.4 


0.18 


WN 


2 


20.77 


20.79 


-4.2 


-4.3 


Fid 


WRIOO 


MJ B4 


J013352.83+304347.8 


0.15 


WNE 


2 


20.56 


20.41 


-4.4 


-4.6 


65 


WRIOI 


MC 58 


J013353.25+304413.6 


0.16 


wc 


2 


20.13 


20.34 


-4.8 


-4.7 


Fid 


WR102 


MC 59 


J013353.60+303851.6 


0.05 


Ofpe/WN9 


2 


18.09 


17.87 


-6.9 


-7.2 


142 


WR103 


MJ X15 



Table 5 — Continued 



Star Type Ref^ mxAi^o^ My^ Mx4750^ OB^ Mjs ID^ Comment 



J013353 


80+303528 

tKJ\J \ KJ\J*JKJ i-J\J > 


,7 


0. 


,17 


WNE 


2 


19 


.41 


18, 


.76 


-5. 


,6 


-6, 


.3 


3 


WR104 


M.1 G34 

±.\3.kJ KJ^^ 


J013354 


40+303453 


,2 


0. 


.20 


WC6 


6 


21, 


.61 


22, 


.35 


-3. 


,4 


-2, 


.7 


Fid 


WR105 


AM 14 


J013354 


85+303222 


.8 


0. 


.29 


WN9-10 


6 


18, 


.34 


18, 


.38 


-6. 


,6 


-6, 


.7 


6 


WR106 

¥ ¥ -I. ly-L WW 


OB 6-5 


.1013355 


33+302001 


.0 


0. 


.76 


WN3+neb 


1 


20, 


.69 


20, 


.32 


-4. 


,3 


-4, 


.7 


111 






J013355 


60+304534 


,9 


0. 


,21 


WN 


2 


21 


.62 


20, 


.92 


-3. 


,4 


-4, 


.1 


67 


WR107 


MJ B17 


,7013355 


71+304501 

• 1 -L I KJKJ^^^KJJ-- 


,3 


0. 


,19 


WN7 


2 


18, 


.60 


18, 


.41 


-6. 


,4 


-6, 


.6 


67 


WR108 

¥ ¥ -L li-L WLJ 


M.1 B8 

d.\A.Kj J—^\J 


,7013355 


87+304528 

i\J 1 1 KJ\J^.KJ I 


,4 


0. 


,20 


WN7 


4 


17, 


.71 


18, 


.41 


-7, 


,3 


-6, 


.6 


67 






.1013355 


94+302732 

tiJ^C \ KJ\J ^ 1 KJ^i 


1 


0. 


,48 


WN3+B2I 


1 


20, 


.61 


18, 


.41 


-4. 


,4 


-6, 


.6 


107 






.1013355 


94+303407 


,6 


0. 


,24 


WC4-5 


2 


19, 


.14 


18, 


.43 


-5. 


,8 


-6. 


.6 


4 


WR109 

¥ » J. tj-L Wty 


AM 15 


,7013356 


23+303241 

1 ^ KJ 1 KJ yj KJ ^ -L- 1 


,6 


0. 


,29 


wc 


2 


20, 


.99 


21, 


.50 


-4. 


,0 


-3, 


,5 


5 


WRllO 


AM 16 


,7013356 

kJ \J J-\JKJ\J\J ■ 


33+303420 

iKJKJ 1 KJ\JKJ^.i-l\J • 


,8 


0. 


,23 


WN3 


1 


20, 


.88 


21, 


.20 


-4. 


,1 


-3, 


,8 


4 






,7013356 

*J \J J~KJKJ'-^\J ■ 


38+303455 

iKJ*^ 1 KJKJKJ^^^J I 


,5 


0. 


,21 


WC5 


6 


20, 


.83 


21, 


.05 


-4. 


,1 


-4, 


,0 


4 


WRlll 


AM 17 


,7013357 

fJ\J^KJKJ\J 1 1 


20+303512 

1 ^\J 1 KJ\JKJKJJ-^< 


,0 


0. 


,21 


WNL 


2 


20, 


.39 


20, 


.28 


-4. 


,6 


-4, 


,8 


107 


WR112 


MCA 14 


.1013357 

%J \J J-KJ\J\J t • 


75+303416 

• 1 KJ \ KJ\JKJ^J-\Ji 


,3 


0. 


,24 


WN 


2 


21, 


.49 


21, 


.89 


-3. 


,5 


-3. 


.2 


4 


WR113 


MC 63 


JUio*jOo. 




A 

,^ 


n 
u. 


OA 


VV IN i-i 




on 


.o4 


ZU, 


. io 


A 


a 
,o 




Q 


A 


VV JTLl 1^ 




J013358, 


.69+303526. 


,5 


0. 


,21 


Bl la+WNE 


2 


16, 


.85 


17, 


.05 


-8. 


,1 


-8, 


.0 


2 


WR115 


OB 2-4 


J013359, 


.39+303337. 


,5 


0. 


,28 


WC6 


2 


19, 


.82 


19, 


.86 


-5. 


,2 


-5, 


.2 


5 


WR116 


AM 19 


J013359, 


.60+303435. 


,4 


0. 


,25 


WN5 


1 


19, 


.39 


19, 


.27 


-5. 


,6 


-5, 


,8 


4 






J013359, 


.79+305150. 


,0 


0. 


,43 


WN2+09I 


1 


18, 


.87 


19, 


.27 


-6. 


,1 


-5, 


,8 


36 


WR153 


UIT274 


J013359, 


.81+303407. 


,2 


0. 


,27 


WNL 


2 


20, 


.41 


20, 


.21 


-4. 


,6 


-4. 


.8 


4 


WR117 


AM 20 


J013400, 


.57+303809. 


,3 


0. 


,15 


WN? 


6 


20, 


.04 


20, 


.04 


-4. 


,9 


-5. 


.0 


Fid 


WR118 


MJ X16 


J013400, 


.90+303918. 


,2 


0. 


,12 


WN5 


2 


20, 


.33 


20, 


.14 


-4. 


,6 


-4, 


,9 


Fid 


WRllO 


MJ X8 


J013401, 


.30+304004. 


,7 


0. 


,12 


WC6 


6 


18, 


.84 


18, 


.87 


-6. 


,1 


-6, 


,2 


Fid 


WR120 


MJ X19 


J013401, 


.73+303620. 


,1 


0. 


,21 


WC4+abs 


6 


19, 


.09 


19, 


.44 


-5. 


,9 


-5. 


.6 


1 


WR121 


MJ X9 


J013402, 


.30+303749. 


,8 


0. 


,17 


WN3-4+abs 


6 


19, 


.42 


19, 


.17 


-5. 


,6 


-5. 


.9 


Fid 


WR122 


UIT 289 



Table 5 — Continued 



Star 




Type 


Ref*' 








Ma4750' 


OB^ 




ID** Comment 


J013402.93+305126.2 


0.41 


WN3+neb 


1 


19.25 


19.01 


-5.7 


-6.0 


36 






J013406.80+304727.0 


0.28 


WN7+neb 


1 


17.20 


17.49 


-7.8 


-7.6 


71 


WR123 


UIT 303 


J013407.85+304145.1 


0.18 


WN 


2 


20.96 


20.80 


-4.0 


-4.2 


94 


WR124 


MC 67 


J013408.23+305234.2 


0.44 


WN4 


1 


20.28 


20.75 


-4.7 


-4.3 


37 






J013408.90+304732.0 


0.29 


WN3 


1 


21.55 


21.66 


-3.4 


-3.4 


71 






J013409.12+303907.2 


0.22 


wc 


2 


20.38 


20.73 


-4.6 


-4.3 


96 


WR125 


MC 68 


J013410.72+305240.7 


0.45 


WN3 


1 


21.30 


21.21 


-3.7 


-3.8 


39 






J01341 1.14+304637.5 


0.27 


WC4 


1 


20.75 


21.32 


-4.2 


-3.7 


75 






J013415.38+303423.2 


0.41 


WN 


2 


20.26 


19.95 


-4.7 


-5.1 


101 


WR126 


AM 24 


J013415.51+303345.0 


0.43 


WN(E) 


2 


19.92 


19.95 


-5.1 


-5.1 


101 


WR127 


AM 26 


J013415.73+303400.7 


0.42 


WNL 


2 


18.67 


18.56 


-6.3 


-6.5 


101 


WR128 


AM 27 


J013415.85+305522.9 


0.54 


WN4 


1 


20.75 


20.67 


-4.2 


-4.4 


42 






J013416.07+303642.1 


0.36 


Ofpe/WN9 


1 


17.95 


17.89 


-7.0 


-7.2 


100 






J013416.28+303646.4 


0.36 


WC6+neb 


1 


19.94 


20.07 


-5.0 


-5.0 


100 


WR129 


MC 70 


J013416.35+303712.3 


0.35 


WN7 


2 


18.13 


17.93 


-6.8 


-7.1 


100 


WR130 


UIT 343 


J013417.17+303250.8 


0.47 


WC4-5 


2 


20.57 


21.01 


-4.4 


-4.0 


101 


WR131 


MC 71 


J013417.21+303334.7 


0.45 


WN3+neb 


1 


20.27 


20.43 


-4.7 


-4.6 


101 






J013418.37+303837.0 


0.34 


Ofpe/WN9 


1 


18.22 


18.02 


-6.8 


-7.0 


(100) 






J013418.74+303411.8 


0.45 


Ofpe/WN9 


2 


19.58 


19.33 


-5.4 


-5.7 


101 


WR132 


UIT 349 


J013419.16+303127.7 


0.54 


WN3 


1 


21.26 


21.09 


-3.7 


-4.0 


Fid 






J013419.58+303801.5 


0.37 


WN7 


1 


20.17 


20.03 


-4.8 


-5.0 


100 






J013419.68+303343.0 


0.48 


WN3-4+ne 


1 


20.67 


20.52 


-4.3 


-4.5 


101 






J013421. 21+303758.2 


0.39 


WN5+neb 


1 


21.12 


21.15 


-3.9 


-3.9 


(100) 






J013421.97+303314.3 


0.51 


WN4 


1 


20.68 


19.97 


-4.3 


-5.1 


102 






J013422.37+303313.7 


0.52 


WN3+abs 


1 


18.11 


17.82 


-6.9 


-7.2 


102 




KlO-03 



Table 5 — Continued 



Star 




Type 


Ref*' 










OB^ 




ID'' 


Comment 


J013422. 54+303317.1 


0.52 


WN 


2 


20.69 


20.49 


-4.3 


-4.6 


102 


WR133 


MC 72 




J013423.02+304650.0 


0.37 


WN4 


1 


19.43 


19.75 


-5.5 


-5.3 


Fid 








J013425.11+301950.3 


0.98 


WN3 


1 


21.01 


20.70 


-4.0 


-4.3 


Fid 








J013426.96+305256.7 


0.51 


WN4.5 


1 


19.75 


19.53 


-5.2 


-5.5 


40 








J013427.30+305229.2 


0.50 


WN3+neb 


1 


21.32 


21.24 


-3.7 


-3.8 


40 








J013429.56+304145.0 


0.43 


WC4 


1 


20.95 


21.39 


-4.0 


-3.7 


(92) 








J013431.45+305716.7 


0.64 


WN3+neb 


1 


20.10 


20.08 


-4.9 


-5.0 


44 


WR134 


MC 73 




J013432. 11+304705.8 


0.45 


WN6 


5 


18.16 


18.03 


-6.8 


-7.0 


84 




N604-WR8 




J013432.18+304903.7 


0.47 


WC4 


1 


20.73 


21.57 


-4.2 


-3.5 


Fid 








J013432. 24+304702. 7 


0.45 


Of/WNL 


5 


17.27 


21.57 


-7.7 


-3.5 


84 




N604-V1 




J013432.32+304659.7 


0.45 


WCE 


2 


17.47 


21.57 


-7.5 


-3.5 


84 


WR135 


MC 74 


N604-WR1 


J013432.50+304703.5 


0.46 


WNIO 


5 


17.86 


21.57 


-7.1 


-3.5 


84 




N604-WR12 




J013432.60+304705.9 


0.46 


WNL 


2 


16.47 


21.57 


-8.5 


-3.5 


84 


WR136 


MC 75 


N604-WR4 


J013432.76+304717.2 


0.46 


Ofpe/WN9 


4 


19.09 


21.57 


-5.9 


-3.5 


84 








J013432.78+304703.5 


0.46 


WC6 


5 


19.13 


21.57 


-5.8 


-3.5 


84 




N604-WR5 




J013433.22+310019.3 


0.73 


WN7+neb 


1 


18.16 


18.42 


-6.8 


-6.6 


45 








J013433.62+304704.6 


0.47 


WN8+neb 


2 


16.83 


18.47 


-8.1 


-6.6 


84 


WR137 


MC 76 


N604-WR6 (blend) 


J013433.82+304656.3 


0.47 


WC4 


5 


19.54 


19.70 


-5.4 


-5.3 


84 




N604-WR7 




J013434.07+304654.3 


0.47 


WNE 


1 


21.52 


19.70 


-3.5 


-5.3 


84 




N604-WR11 




J013438.18+304953.2 


0.53 


WN5 


1 


20.92 


20.89 


-4.1 


-4.2 


Fid 








J013438.98+304119.8 


0.54 


WC4 


1 


20.20 


20.60 


-4.8 


-4.4 


90 








J013440.42+304321.9 


0.54 


WN4+neb 


1 


19.24 


19.05 


-5.7 


-6.0 


87 


WR138 


MCA 16 




J013443.51+304919.4 


0.58 


WN4 


1 


20.88 


20.62 


-4.1 


-4.4 


Fid 








J013444.28+303757.2 


0.66 


WN3 


1 


21.40 


21.24 


-3.6 


-3.8 


98 








J013444.61+304445.4 


0.58 


WC5-6 


1 


20.31 


20.50 


-4.7 


-4.5 


86 


WR139 


MC 77 





Table 5 — Continued 



Star 




Type 


Ref^ 










OB^ 




ID^ Comment 


J013447.32+310748.8 


1.00 


WN3 


1 


20.79 


20.98 


-4.2 


-4.1 


Fid 






J013458.89+304129.0 


0.78 


WC4 


2 


20.66 


21.28 


-4.3 


-3.8 


88 


WR140 


MC 78 


J013505.37+304114.9 


0.86 


WC4-5 


2 


19.06 


19.12 


-5.9 


-5.9 


89 


WR141 


MC 79 


J013507.24+304500.9 


0.84 


WC5 


1 


20.30 


21.14 


-4.7 


-3.9 


Fid 






J013509.73+304157.3 


0.90 


Ofpe/WN9 


1 


18.04 


17.94 


-6.9 


-7.1 


89 




Romano's star 


J013510.27+304522.9 


0.87 


WN2pec 


1 


21.81 


22.02 


-3.2 


-3.0 


Fid 







^Distance from the center within the plane of M33, normahzed to the D-2b isophotoal radius of 30f8 assuming 0:2000 = 
0l'*33"*50?89, ^2000 = 30°39'36'.'8, an inclination of 56° and a position angle of the major axis of 23°, following Kwitter & 
Aller 1981 and Zaritsky et al. 1989. 

'^References for spectral types: 1 — This work; 2 — Massey & Johnson 1998; 3 — Kehrig et al. 2011; 4 — Massey et al. 2007b; 
5— Drissen et al. 2008; 6— Abbott et al. 2004 

*^From Massey et al. 2005. 

^AB magnitude through CT filter centered at 4750 A and calibrated using the values from Massey & Johnson 1998. 

*^ Absolute magnitudes computed assuming a true distance modulus of 24.60 (830 Mpc) and adopting an average reddening 
of E[B — V) = 0.12 based on Massey et al. 2007a. Adopting Ry = 3.1 leads to an Ay = 0.37 mag and A^^^^^ = 0.45 mag. 

^OB association as defined by Humphreys & Sandage 1980. Parenthesis means the star is just outside the boundaries of 
the association, while "Fid" implies it is a field star, not in a cataloged OB association. 

s Massey &; Johnson 1998 WR designation. 

''IDs are from Massey k, Johnson 1998 and references therein, other than the HII regions designations N592 and N604 



from Drissen et al. 2008 and the "KIO" designations from Kehrig et al. 2011). 

^J013332.64+304127.2 was shown to be a late-type WN star by Massey & Conti 1983 (MC28), and later identified by 
Hell A4686 imaging by Armandroff k, Massey 1985, where it is listed as AM 2. However, Massey et al. 2006 found that the 
spectrum has since changed dramatically, resembling a hot LBV. 
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Table 6. Isolated WR Stars in M33 



Star 


P 


Type 


J013425.11+301950.3 


0.98 


WN3 


J013438.18+304953.2 


0.53 


WN5 


J013507.24+304500.9 


0.84 


WC5 


J013510.27+304522.9 


0.87 


WN2pec 



Table 7. Absolute Visual Magnitudes'' 



Subtype 




M33 






LMC 






SMC 




My 


u 




My 


a 


N 


My 


a 


N 


WN2 


-6.0 


0.9 


3 


-2.6 




1 


-6.8 




1 


WN3 


-4.8 


1.2 


25 


-4.0 


0.6 


11 


-5.0 


1.2 


6 


WN4-4.5 


-5.2 


1.3 


19 


-4.2 


0.6 


33 


-5.2 




1 


WN5 


-5.4 


1.5 


5 


-5.1 


1.4 


3 








WN6 


-4.9 


0.9 


5 


-4.8 


1.5 


3 


-6.2 




1 


WN7 


-5.6 


1.4 


10 


-5.9 


1.0 


4 








WN8 


-4.5 


0.9 


6 


-5.5 


0.0 


2 








WN9 


-4.4 


0.9 


11 


-5.2 


0.7 


2 








WC4 


-5.9 


1.3 


14 


-5.0 


0.9 


11 








WC5 


-4.9 


0.8 


4 














WC6 


-4.5 


0.9 


11 














WC7 


-4.2 


0.3 


2 















''The values are the means of the WRs of a given subtype excluding ones 
with uncertain types and/or with companions. 



-61 - 



Table 8. Additional WR Candidates 



Star 


V 


"^A4650 


J013302.73+301131.6 


20.00 


19.65 


J013404.07+304658.3 


19.53 


19.40 


J013411.45+303637.3 


21.34 


21.24 


J013432.60+304211.3 


19.34 


19.10 


J013434.26+304637.8 


20.48 


20.71 


J013442.41+305019.0 


20.95 


21.19 



Table 9. Numbers of WC and WN Stars 



Region 


P 


log(0/H) + 12 


#WCs 


#WNs 


WC/WN 


p < 0.25 


0.16 


8.72 


26 


45 


0.58±0.09 


0.25 >p< 0.50 


0.38 


8.41 


15 


54 


0.28±0.07 


p > 0.50 


0.69 


8.29 


12 


54 


0.22±0.06 


LMC 




8.37 


25 


109 


0.23±0.06 


SMC 




8.13 


1 


11 


0.09±0.05 



